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Abstract 



The energy and time scales during the inflationary stage of the universe calls 
for an out of equilibrium quantum field treatment. Moreover, the high energy 
■ densities involved (~ 1/g ~ 10^^) make necessary the use of non-perturbative 

ly^ . approaches as the large N and Hartree methods. We start these lectures by 

'nI" I introducing the such non-perturbative out of equilibrium methods in cosmo- 

I logical universes. We discuss the renormalization procedure and the choice 

OO ' of initial conditions. We then study with these methods the non-linear dy- 

namics of quantum fields in matter and radiation dominated FRW and de 
Sitter universes. For a variety of initial conditions, we compute the evolu- 
, tion of the inflaton, its quantum fluctuations and the equation of state. We 

Q I investigate the phenomenon of explosive particle production due to spinodal 

' instabilities and parametric amplification in FRW and de Sitter universes 

with and without symmetry breaking. We find that the particle production 
^ . is somewhat sensitive to the expansion of the universe. In the large N limit 

for symmetry breaking scenarios, we determine generic late time fields be- 
havior for any flat FRW and de Sitter cosmology. We find that quantum 
fluctuations damp in FRW as the square of the scale factor while the or- 
der parameter approaches a minimum of the potential at the same rate. We 
present a complete and numerically accessible renormalization scheme for the 
equation of motion and the energy momentum tensor in flat cosmologies. In 
this scheme the renormalization constants are independent of time and of the 
initial conditions. Furthermore, we consider an 0{N) inflaton model coupled 
self-consistently to gravity in the semiclassical approximation, where the field 
is subject to 'new inflation' type initial conditions. We study the dynamics 
self-consistently and non-perturbatively with non-equilibrium field the- 
ory methods in the large N limit. We find that spinodal instabilities drive 
the growth of non-perturbatively large quantum fluctuations which shut off 
the inflationary growth of the scale factor. We find that a very specific com- 
bination of these large fluctuations plus the inflaton zero mode assemble into 
a new effective field. This new field behaves classically and it is the object 
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which actually rolls down. We show how this reinterpretation qualitatively 
saves the standard picture of how metric perturbations are generated dur- 
ing inflation and that the spinodal growth of fluctuations dominates the time 
dependence of the Bardeen variable for superhorizon modes during inflation. 
We compute the amplitude and index for the spectrum of scalar density and 
tensor perturbations and argue that in all models of this type the spinodal 
instabilities are responsible for a 'red' spectrum of primordial scalar density 
perturbations. The decoherence aspects and the quantum to classical transi- 
tion through inflation are studied in detail by following the full evolution of 
the density matrix. 



I. INTRODUCTION AND MOTIVATION 

Inflationary cosmology has come of age. From its beginnings as a solution to deep 
problems such as the horizon, flatness, entropy and monopole problems it has grown 
into the main contender for the source of primordial fluctuations giving rise to large scale 



structure ||33|. There is evidence from the measurements of temperature anisotropics in the 
cosmic microwave background radiation (CMBR) that the scale invariant power spectrum 
predicted by generic inflationary models is at least consistent with observations [0-^ and we 
can expect further and more exacting tests of the inflationary power spectrum when the MAP 
and PLANCK missions are flown. In particular, if the fluctuations that are responsible for 
the temperature anisotropics of the CMB truly originate from quantum fluctuations during 
inflation, determinations of the spectrum of scalar and tensor perturbations will constrain 
inflationary models based on particle physics scenarios and probably will validate or rule 
out speciflc proposals 

The tasks for inflationary universe researchers are then two-fold. First, models of inflation 
must be constructed on the basis of a realistic particle physics model. This is in contrast 
to the current situation where most, if not all acceptable inflationary models are ad- hoc in 
nature, with flelds and potentials put in for the sole purpose of generating an inflationary 
epoch. Second, and equally important, we must be sure that the quantum dynamics of 
inflation is well understood. This is extremely important, especially in light of the fact 
that it is exactly this quantum behavior that is supposed to give rise to the primordial 
metric perturbations which presumably have imprinted themselves in the CMBR. This latter 
problem is the focus of this review. 

The inflaton must be treated as a n on- equilibrium quantum fleld . The simplest way to 
see this comes from the requirement of having small enough metric perturbation amplitudes 
which in turn requires that the quartic self coupling A of the inflaton be extremely small, 
typically of order ~ 10"^^. Such a small coupling cannot establish local thermodynamic 
equilibrium (LTE) for all fleld modes; typically the long wavelength modes will respond 
too slowly to be able to enter LTE. In fact, the superhorizon sized modes will be out of 
the region of causal contact and cannot thermalize. We see then that if we want to gain 
a deeper understanding of inflation, non-equilibrium tools must be developed. Such tools 
exist and have now been developed to the point that they can give quantitative answers to 



these questions in cosmology ^ - [0, P^j27|j3^j40[] . These methods permit us to follow the 



dynamics of quantum flelds in situations where the energy density is non-perturbatively 
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large (~ 1/A). That is, they allow the computation of the time evolution of non-stationary 
states and of non-thermal density matrices. 

Our programme on non-equilibrium dynamics of quantum field theory, started in 1992 
0, is naturally poised to provide a framework to study these problems. The larger goal 
of the program is to study the dynamics of non-equilibrium processes from a fundamental 
field-theoretical description, by solving the dynamical equations of motion of the underlying 
four dimensional quantum field theory for physically relevant problems: the early universe 
dynamics, high energy particle collisions, phase transitions out of equilibrium, symmetry 
breaking and dissipative processes. 

The focus of our work is to describe the quantum field dynamics when the energy density 
is high. That is, a large number of particles per volume m~^, where m is the typical mass 
scale in the theory. Usual S-matrix calculations apply in the opposite limit of low energy 
density and since they only provide information on in out matrix elements, are unsuitable 
for calculations of expectation values. 

In high energy density situations such as in the early universe, the particle propagator (or 
Green function) depends on the particle distribution in momenta in a nontrivial way. This 
makes the quantum dynamics intrinsically nonlinear and calls to the use of self-consistent 
non-perturbative approaches as the large N limit, Hartree and self-consistent one- loop ap- 
proximations. 

There are basically three different levels to study the early universe dynamics: 

1. To work out the nonlinear dynamics of quantum fields in Minkowski spacetime. By 
non-linear dynamics we understand to solve the quantum equations of motion includ- 
ing the quantum back-reaction quantitatively - [0, p4| , |3ll , ^^ . This level is in 



fact appropriate to describe high energy particle collisions |^ . 
2. To work out the nonlinear dynamics of quantum fields in fixed cosmological back- 



grounds 1 14, 15]. New phenomena arise then compared with 1. showing that a 



Minkowski analysis is not quantitatively precise for expanding universes. 

3. A self-consistent treatment of the quantum fields and the cosmological background 
[0,0. That is, the metric is obtained dynamically from the quantum fields (matter 
source) propagating in the that metric. 

We shall successively present the three levels of study. The first stage was reviewed 
in the 1996 Chalonge School |]TT| (see []T2[ for further progress). The second level is the 



subject of sees. Vll and VIII. We study the parametric and spinodal resonances both in 
FRW and de Sitter backgrounds wide range of initial conditions both in FRW and de Sitter 
backgrounds |13,1^. [Parametric resonance appears in chaotic inflationary scenarios for 



unbroken symmetry whereas spinodal unstabilities show up in new inflation scenarios with 
broken symmetry]. Both types of unstabilities shut-off through the non- linear quantum 



evolution as described in sees. VII and VIII [13,15] both analytically and numerically. We 



follow the equation of state of the quantum matter during the evolution and analyze its 
properties. 

The third stage of our approach is to apply non-equilibrium quantum field theory tech- 
niques to the situation of a scalar field coupled to semiclassical gravity, where the source of 
the gravitational field is the expectation value of the stress energy tensor in the relevant. 



3 



dynamically changing, quantum state. In this way we can go beyond the standard analyses 
||35|^8| which treat the background as fixed and do not consider the non-linear quantum 
field dynamics. 

In all cases 1. - 3. , the quantum fields energy- momentum tensor is covariantly conserved 



both at the regularized as well as the renormalized levels [|I^ - [|T^. 

We mainly consider for the stage 3. new infiation scenarios where a scalar field <^ evolves 
under the action of a typical symmetry breaking potential. The initial conditions will be 
taken so that the initial value of the order parameter is near the top of the potential (the 
disordered state) with essentially zero time derivative. What we find is that the existence 
of spinodal instabilities, i.e. the fact that eventually (in an expanding universe) all modes 
will act as if they have a negative mass squared, drives the quantum fiuctuations to grow 
non-perturbatively large. We have the picture of an initial wave-function or density matrix 
peaked near the unstable state and then spreading until it samples the stable vacua. Since 
these vacua are non-perturbatively far from the initial state (typically ~ m/V^, where m 
is the mass scale of the field and A the quartic self-coupling), the spinodal instabilities will 
persist until the quantum fluctuations, as encoded in the equal time two-point function 
($(f, t)^), grow to 0{im?/X). 

This growth eventually shuts off the inflationary behavior of the scale factor as well as the 
growth of the quantum fluctuations (this last also happens in Minkowski spacetime [10-0). 

The scenario envisaged here is that of a quenched or super-cooled phase transition where 
the order parameter is zero or very small. Therefore one is led to ask: 

a) What is rolling down?. 

b) Since the quantum fluctuations are non-perturbatively large ( ~ 1/A), will not they 
modify drastically the FRW dynamics?. 

c) How can one extract (small?) metric perturbations from non-perturbatively large fleld 
fluctuations? 

We address the questions a)-c) as well as other issues in sec. IX. 

We choose such type of new inflationary scenario because the issue of large quantum fluc- 
tuations is particularly dramatic there. However, our methods do apply to any inflationary 
scenario as chaotic, extended and hybrid inflation. 



II. QUANTUM FIELD THEORY AROUND AN EXCITED STATE WITH 
NON-PERTURBATIVE ENERGY DENSITY 

To start we present a simple study of quantum fleld theory around an excited state. 
The relevant situation both in cosmology and high energy particle scattering is when the 
energy density is large (~ C(l/A), non-perturbative) with respect to the ground state. 
For simplicity we hall make the derivation in Minkowski spacetime. The generalization to 
cosmological spacetimes is straightforward. 

Let us consider a scalar fleld model with quartic selfcoupling 

^ = \{d,^f-'^^'-j,^'- (2.1) 
Here $(a;) is a real fleld. The equations of motion take the form 
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9^$ + $ + - $ 
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and the canonical momentum is 



U(x) 



(2.2) 



(2.3) 



For classical fields, the solutions of the equations of motion ( p. 21) depend on the coupling 
constant A as follows, 



^c{'x) = —j= Fimx) 
V A 

where the index ^ indicates that we are considering c-number solutions of eqs. ( p.2| ) and F is 
a dimensionless functions. Therefore, the classical solutions are large for small coupling and 
the energy density 



«=2 



1 r- 



+ (V$)^ 



H $^ + — <l>^ 

2 4! 



scales as m^/X and is large too. Notice that here X ^ h. Therefore, for small A we should 
expect some semiclassical behavior. 

Let us now consider the quantum field $(x). $(x) and n(a;) will be operators in a Fock 
space obeying the canonical commutation rules. 



[<i>(f,t),n(y,t)] =^5(f-y) 



(2.4) 



Let us consider a quantum state | > which is not the vacuum. The expectation value of 
^(x) there 

0(x) = < <l>(x) > 

will in general be a function of space and time, (pix) will be typically of order for small 
A just for correspondence with the classical theory. 
We write the operator $(a;) as, 

$(x) = 0(x) + Ipix) , 

where iplx) is a new quantum operator. iIj{x) must obey the constraint 

< ipix) >= . 



The equations of motion ( |2.2| ) hold in the quantum theory for the operator $(x). Taking 
the expectation value of eq. (|2.2| ) in the state | > yields 



(2.5) 



((9^ + m^) (f){x) + ^(l)ixf + ^(f){x) < ^pix)'^ > +7: < >= 

^ ^ o 2 



Subtracting eq. (|2.5|) from eq.(p.2|) yields. 



ijj{x) + ^(f){x) V'(x)^- < V^(x)^ > 



+ 



A 



ipix)^- < ijj{xf > 



. 

(2.6) 
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In order to solve (|2.6| ) perturbatively in A but keeping (f){x) arbitrary, we introduce the 
Green function, 



9^ + + ^(f){xf 



G{x,y) = S^{x - y) . 



We consider the retarded Green function which vanishes for x^ < y^. This Green function 
(propagator) depends on the field expectation value (j){x) showing that the particle propa- 
gation depends here on the properties of the state. As is well known, this is not the case on 
the vacuum where the propagator is explicitly known. 

With the help of this retarded Green function, we can write the exact evolution equation 
( p.6|) as follows 

^{x) = Mx) G(x, y) {30(y) [i;{yf- < i;{yf >] + i;{yf- < ij{yf >} (2.7) 

Here, il)o{x) stands for the general operator solution of the homogeneous equation: 



i)Q{x) = 



with < tpo{x) >= 0. ipoix) will be expressed as a sum of c-number eigenfunctions times 
creation and annihilation operators (see below eq. ( |2.10| ). It then follows that 

< ipo{x)'^ >^ but < ipoixf >= . 

For small A we see from eq. (|2.7|) that 

ij{x) = M^) + o{X) . 

Inserting that information into eq. ( |2.5| ) yields, 



9^ + ^ < ipoixY > 



0o(x) + XM^Y = 



(2.8) 



where, 



0(x) = 0o(x) + 0(A) . 
In addition, we can consider ipo^x) as the general operator solution of the equation 



(2.9) 



Let us consider from now on spatially homogeneous states | >. Then, 

< $(x) >= (j){t) 

only depends on time. This is the relevant situation for early universe investigations. We 
can now Fourier expand ipo{x) as follows. 



ipo{x,t) 



v^(27r)3/2 



(2.10) 
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Here, and stand for creation and annihilation operators on the state | >. That is. 



agl >= and at| >^ for all k 



and 



(2.11) 



(2.12) 



The mode functions /^(t) have here a non-trivial time dependence. [Recall that they have 
just harmonic time dependence around the vacuum]. Actually, most of the field dynamics 
is encoded in their evolution defined by eq. (p.9|) 



fk(t) = 

The canonical momentum ( |2.3| ) is Fourier expanded as follows, 

d^k 



TTo{x, t) 



v^(27r)3/2 



-ik-x 



(2.13) 



It is easy to show now that the canonical commutation rules ( p^.4D hold using eg. ( p^. 121 ) and 
the time independence of the Wronskian that we normalize as follows: 

It is then convenient to choose as initial conditions for the mode functions: 

1 



fk{0) 



fk{0) 



(2.14) 



with 



'A;2 + m2 + ^ 0(0)2 



One must specify in addition the values of 

0(0) = 00 and 0(0) = po . 



(2.15) 



We can now compute the expectation value < ipo{x)^ > that appears in eq.( |2.8| ). We get 
using eqs. (|2l0|) , (|Tl|) and (|T|, 



< >-- 



d^k 

(2^ 



P dk 
27r2 



This integral diverges in the ultraviolet. We put for the moment a cutoff A in the momentum. 
We shall discuss below (sec. V) how the cutoff dependence can be absorbed into mass a 
coupling constant renormalization. 

In summary, we have obtained a set of self-consistent evolution equations for the quantum 

state I >: 
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fk{t) = 



27r2 



\fk'{t)\' 



(2.16) 



with the initial conditions defined by eqs.( 2.14 )-( 2.15 ). The choice of the initial conditions 
defines the state | >. We have an infinite number of unknowns: 

/fc(t),0o(^) for < A; < oo and t>0. 



defined by the coupled nonlinear equations ( |2.16D . 

All physical quantities can be computed in terms of the mode functions fk(t), < k < oo 
and the order parameter 0o(^)- We have described a small coupling or self-consistent one- 
loop approximation. The large N and Hartree approximations are discussed in sec. IV. 
Large N and Hartree contain definitely more information than the self-consistent one-loop 
approximation (|2.16|) . However, the simplicity of its derivation makes it a very useful and 
pedagogical exercise. 



III. NON-EQUILIBRIUM QUANTUM FIELD THEORY, SEMICLASSICAL 

GRAVITY AND INFLATION 

We present here the framework of the non-equilibrium closed time path formalism. For 
a more complete discussion, the reader is referred to |T^, or the alternative approach given 
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The time evolution of a system is determined in the Schrodinger picture by the functional 
Liouville equation 

'^ = iH{t),p{t)], (3.1) 

where p is the density matrix and we allow for an explicitly time dependent Hamiltonian as is 
necessary to treat quantum fields in a time dependent background. Formally, the solutions to 
this equation for the time evolving density matrix are given by the time evolution operator, 
U{t, t), in the form 

p{t) = U{t,to)p{to)U-\t,to). (3.2) 

The quantity p(to) determines the initial condition for the evolution. We choose this initial 
condition to describe a state of local equilibrium in conformal time, which is also identified 
with the conformal adiabatic vacuum for short wavelengths. In the appendix we provide an 
analysis and discussion of different initial conditions and their physical content within the 
context of expanding cosmologies. 

Given the evolution of the density matrix (|3.2|) , ensemble averages of operators are given 
by the expression (again in the Schrodinger picture) 

^ Tr[Uito,t)OUit,t')Uit',to)p{to)] ^ 
Trp{to) 
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where we have inserted the identity, U{t,t')U{t' ,t) with t' an arbitrary time which will be 
taken to infinity. The state is first evolved forward from the initial time to to t when the 
operator is inserted. We then evolve this state forward to time t and back again to the 
initial time. 

The actual evolution of various quantities in the theory can now be evaluated by either 
constructing the appropriate Green functions as in [jTO[, or by choosing an explicit Ansatz 
for the functional form of the time dependent density matrix so that the trace in (|3.3| ) may 
be explicitly evaluated as a functional integral (see [0). The methods are equivalent, and 
provide the results which will be presented below for the cases of interest. 

We shall study the infiationary dynamics in a spatially flat Friedmann-Robertson- Walker 
background with scale factor a{t) and line element: 

ds'^ = dt^ - a^{t) dx^. (3.4) 

Our Lagrangian density has the form 

1. 



(3.5) 



Our approach can be generalized to open as well as closed cosmologies. 

Our program incorporates the non-equilibrium behavior of the quantum fields involved 
in infiation into a framework where the geometry (gravity) is dynamical and is treated self 
consistently. We do this via the use of semiclassical gravity where we say that the 
metric is classical and determined through the Einstein equations using the expectation 
value of the stress energy tensor {T^y). Such expectation value is taken in the dynamically 
determined state described by the density matrix p{t). This dynamical problem can be 
described schematically as follows: 

1. The dynamics of the scale factor a{t) is driven by the semiclassical Einstein equations 

o \ + J^^ g^y + (higher curvature) = -{T^u)r- (3.6) 

Here Gr, are the renormalized values of Newton's constant and the cosmological 
constant, respectively and G^j, is the Einstein tensor. The higher curvature terms 
must be included to absorb ultraviolet divergences. 

2. On the other hand, the density matrix p(t) of the matter (that determines (T^i/)ij) 
obeys the Liouville equation 

^^ = [H,p{t)], (3.7) 

where H is the evolution Hamiltonian, which is dependent on the scale factor, a{t). 

It is this set of equations we must try to solve; it is clear that initial conditions must be 
appended to these equations for us to be able to arrive at unique solutions to them. Let us 
discuss some aspects of the initial state of the field theory first. 
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A. On the initial state: dynamics of phase transitions 



The situations we consider are 

1. the theory admits a symmetry breaking potential and in which the field expectation 
value starts its evolution near the unstable point. 

2. The symmetry is not broken and the field expectation value starts its evolution at a 
finite distance from the absolute minimum. 

There is an issue as to how the field got to have an expectation value near the unstable 
point (typically at $ = 0) as well as an issue concerning the initial state of the non-zero 
momentum modes. 

Since our background is an FRW spacetime, it is spatially homogeneous and we can 
choose our state p(t) to respect this symmetry. Starting from the full quantum field t) 
we can extract a part that has a natural interpretation as the zero momentum, c-number 
part of the field by writing: 

(Pit) = TT[p{tMx,t)] = mx,t))). (3.8) 

The quantity ^(£, t) represents the quantum fiuctuations about the zero mode (j)(t) and 
clearly satisfies (\l/(x, t)) = 0. 

We need to choose a basis to represent the density matrix. A natural choice consistent 
with the translational invariance of our quantum state is that given by the Fourier modes, 
in comoving momentum space, of the quantum fiuctuations \E'(x, t): 

-^{x, t)= j exp(-z k ■ x) Mt)- (3.9) 

In this language we can state our ansatz for the initial condition of the quantum state 
as follows. We take the zero mode (pit = 0) = 0o; 0(^ = 0) = 0, where 0o "will typically be 
very near the origin for broken symmetry and at a finite distance from it in the unbroken 
symmetry case. The initial conditions on the the nonzero modes ipki^ = 0) will be chosen 
such that the initial density matrix p(t = 0) describes a vacuum state (i.e. an initial state 
in local thermal equilibrium at a temperature Tj = 0). There are some subtleties involved 



in this choice. First, as explained in [|I4[, in order for the density matrix to commute 
with the initial Hamiltonian, we must choose the modes to be initially in the conformal 
adiabatic vacuum (these statements will be made more precise below). This choice has the 
added benefit of allowing for time independent renormalization counterterms to be used in 
renormalizing the theory. 

We are making the assumption of an initial vacuum state in order to be able to proceed 
with the calculation. It would be interesting to understand what forms of the density matrix 
can be used for other more general initial conditions. 

The assumptions of an initial equilibrium vacuum state are essentially the same used in 



refs. [^, and in the analysis of the quantum mechanics of infiation in a fixed de 
Sitter background. 
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As discussed in the introduction, if we start from such an initial state, spinodal or para- 
metric instabihties will drive the growth of non-perturbatively large quantum fluctuations. 
In order to deal with these, we need to be able to perform calculations that take these large 
fluctuations into account. Although the quantitative features of the dynamics will depend on 
the initial state, the qualitative features associated with spinodal or parametric unstabilities 
are fairly robust for a wide choice of initial states that describe a phase transition. 



IV. THE INFLATON MODEL AND THE EQUATIONS OF MOTION 



Having recognized the appearance of large quantum fluctuations driven by parametric or 
spinodal unstabilities, we need to study the dynamics within a non-perturbative framework. 
That is, a framework allowing calculations for non-perturbatively large energy densities. 
We require that such a framework be: i) renormalizable, ii) covariant energy conserving, iii) 
numerically implementable. There are very few schemes that fulflll all of these criteria: the 
Hartree and the large N approximation |T0|- |jl6|. Whereas the Hartree approximation is 



basically a Gaussian variational approximation pSj that in general cannot be consistently 
improved upon, the large N approximation can be consistently implemented beyond leading 
order [^,^]. In addition, the presence of a large number of fields in most of the GUT's 
models suggest that the large limit will be actually a realistic one. Moreover, for the case 
of broken symmetry it has the added bonus of providing many light fields (associated with 
Goldstone modes) that will permit the study of the effects of other fields which are lighter 
than the inflaton on the dynamics. Thus we will study the infiationary dynamics within 
the framework of the large A^ limit of a scalar theory in the vector representation of 0{N) 
both for unbroken and broken symmetry. In the second case we will have a quenched phase 
transition. 



We assume that the universe is spatially flat with a metric given by eq. (|3.4| ). The matter 
action and Lagrangian density are given by eq.(^]5|). 



Sm = j d'^x Cm = j d'^x a^{t) 



•2 



[X] 



1 (V$(x))^ 

2 a2(t) 



- V{^{x)) 



(4.1) 



where m? > for unbroken symmetry and < for broken symmetry. Here 7^(t) stands 
for the scalar curvature 

The coupling of $(a;) to the scalar curvature 7^(t) has been included since arises anyhow as 



a consequence of renormalization |13] 



The gravitational sector includes the usual Einstein term in addition to a higher order 
curvature term and a cosmological constant term which are necessary to renormalize the 
theory. The action for the gravitational sector is therefore: 
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d'^x a^{t) 



167rG 2 ^ ' 



(4.4) 



with K being the cosmological constant (we use K rather than the conventional A/8ttG to 
distinguish the cosmological constant from the ultraviolet cutoff A we introduce to regular- 
ize the theory; see section V). In principle, we also need to include the terms R^" R^i, and 
as they are also terms of fourth order in derivatives of the metric (fourth adia- 
batic order), but the variations resulting from these terms turn out not to be independent 
of that of 7^2 in the flat FRW cosmolo gy we are considering. 

The variation of the action S = Sg + with respect to the metric g^^ gives us Einstein's 
equation 



G 



8nG 



+ aHau + Kg 



-<T.,> , 



(4.5) 



where G^^, is the Einstein tensor given by the variation of y/—gTZ, H^i, is the higher order 
curvature term given by the variation of ^J—gTZ'^, and T^^ is the contribution from the 
matter Lagrangian. With the metric ( p.4| ), the various components of the curvature tensors 
in terms of the scale factor are: 



-?>{d/af , 



' (X- 

-n + —: 

a 



-7^ = -6 


/ d 


- + 








= 


12 / ' 





-6 (7^ + 3-7^ 

a 



Eventually, when we have fully renormalized the theory, we will set = and keep as our 
only contribution to Kji a piece related to the matter fields which we shall incorporate into 
T, 



A. The Large N Approximation 

To obtain the proper large limit, the inflaton field is written as 

$(f,t) = (cT(f,t),7r(f,t)), 
with TT an (iV — l)-plet, and we write 

a{x,t) = VN(j){t) + x{x,t) ; {a{x,t)) = VN<j){t) ; {x{x,t)) = 0. (4.6) 
To implement the large N limit in a consistent manner, one may introduce an auxiliary field 



as m 



the factorization 15,14 



]40[| . However, the leading order contribution can be obtained equivalently by invoking 

Q{X^)X^ + constant , 5{x^)x , 

(tt ■ n f 2{n^)n^ - {nY + Oil/N) , 

^\'^{n')x' + ^'{x'),n\^{n')x. (4.7) 
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To obtain a large limit, we define [p!5|,|14 



N-l 



7l{x,t) = ij{x,t) (1, 1, ■ ■ 1), 

where the large N limit is implemented by the requirement that 



(4.8) 



(4.9) 



The leading contribution is obtained by neglecting the 0{1/N) terms in the N ^ oo limit. 
The resulting Lagrangian density is quadratic, with linear terms in x and vr. The equations of 
motion are obtained by imposing the tadpole conditions < x(a;, t) >= and < 7f(x, t) >= 
which in this case are tantamount to requiring that the linear terms in x ^ iii the 
Lagrangian density vanish. Since the action is quadratic, the quantum fields can be expanded 
in terms of creation and annihilation operators and mode functions that obey the Heisenberg 
equations of motion 



7f(x, t) 



(2^ 



-ik-x 



(4.10) 



We see that since there are A^ — 1 'pion' fields, contributions from the field x can be neglected 
in the iV — >■ cx) limit as they are of order 1/A^ with respect those of ip and 0. 



The tadpole condition leads to the following equations of motion [15,14|: 



^(t) + 3H{t) (/.(t) + M\t) 0(t) = 0, 



with the mode functions 



fc2 



fk{t) = 0, 



where 



\m+\mt)). 



(4.11) 



(4.12) 



(4.13) 



An important point to note in the large A^ equations of motion is that the form of the 
equation for the zero mode (^4.11 ) is the same as for the = mode function ( [4.12 ). It will 
be this identity that allows solutions of these equations in a symmetry broken scenario to 
satisfy Goldstone's theorem. 

In this leading order in 1/A^ the theory becomes Gaussian, but with the self-consistency 
condition 



2(2^ 



(4.14) 



The initial conditions on the modes fk{t) must now be determined. At this stage it proves 
illuminating to pass to conformal time variables in terms of the conformally rescaled fields 
(see and section VI for a discussion) in which the mode functions obey an equation 



which is very similar to that of harmonic oscillators with time dependent frequencies in 
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Minkowski space-time. It has been realized that different initial conditions on the mode 
functions lead to different renormalization counterterms |14]; in particular imposing initial 
conditions in comoving time leads to counterterms that depend on these initial conditions. 
Thus we chose to impose initial conditions in conformal time in terms of the conformally 
rescaled mode functions leading to the following initial conditions in comoving time: 



fkito) 



1 



fkito) 



a{to) 



fk{to), 



with 



W^ = k^ + M\to 



6 



(4.15) 



(4.16) 



For convenience, we have set a(to) = 1 in eq.( [4.16|) . At this point we recognize that when 
-^^(^o) — Tlitoi/Q < the above initial condition must be modified to avoid imaginary 
frequencies, which are the signal of instabilities for long wavelength modes in the broken 
symmetry case. Thus we define the initial frequencies that determine the initial conditions 
( CT ) as 



wi = e + M\to 



7^(to) 



6 
6 



for fc^ < 
for > 



7^(^o) 



7^(^o 



(4.17) 
(4.18) 



In the unbroken symmetry case (m^ > ) we use eq.( 4.18| ) for all k. 

As an alternative we have also used initial conditions which smoothly interpolate from 
positive frequencies for the unstable modes to the adiabatic vacuum initial conditions defined 
by ( [4.15|) -( p:.16|) for the high k modes. While the alternative choices of initial conditions 
result in small quantitative differences in the results (a few percent in quantities which 
depend strongly on these low-fc modes), all of the qualitative features we will examine are 
independent of this choice. 



In the large A^ limit we find the energy density and pressure density to be given by p!5|J14 

,4 



e 



1 > 



— = —m H — m 
A^ 2^ 2 



2 j2 



(4.19) 



e — 3p 
N 



A 

+ 2 



+ 



2m^ 



X 



+ 



+ 3-^ 
a 



(4.20) 



where {ip'^) is given by equation ( |4.14|) and we have defined the following integrals: 
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<Pk 



2(27r)3 



(4.21) 
(4.22) 



The composite operators {ipip) and {ipip) are symmetrized by removing a normal ordering 
constant to yield 



(27r)3 
d^k 
(2^ 



dt 



fk{t)m + ut)m) 



(4.23) 
(4.24) 



The last of these integrals, ( [4.24| ), may be rewritten using the equation of motion ( |4.12| ): 

(4.25) 



a 



It is straightforward to show that the bare energy is covariantly conserved by using the 
equations of motion for the zero mode and the mode functions. 



B. Hartree Approximation 

In the Hartree approximation, our theory is that of a single component scalar field, 
$(x, t), with the Z2 symmetry $ — $. The potential can be written as: 

y($) = i(m2 + e7e)$^ + ^< (4.26) 

where TZ is the Ricci scalar. We have rescaled the coupling constant here by a factor 3 
compared with the one used in the large limit [eg. ( [4.21) ]. 

We decompose the field into its zero mode, (pit) = ($(x, t)), fluctuations ip^x^t) about 

it: 

$(x,t) = 0(t) + V(x,t). 



The potential ( [4.26| ) may then be expanded in terms of these fields. 



The Hartree approximation is achieved by making the potential quadratic in the fluctu- 
ation field ip by invoking the factorization 

^3(x,t) -.3(^2(f,t))^(f,t) , 

This factorization yields a quadratic theory in which the effects of interactions are encoded 
in the time dependent mass which is determined self-consistently. 

The equations of motion for the zero mode and the fluctuations are given by the tadpole 
equation 

(^(f,t)) = 0. 



15 



Introducing the Fourier mode functions, Uk{t), they can be written as: 



a(r) D 2 



cP a{t) d 



Uk{t) = . 



(4.27) 
(4.28) 



2(27r)3 



(4.29) 



The initial conditions on the mode functions are the same as eq. (|4.15|) in the large limit, 



a{to 



Uk{t 



0), 



(4.30) 



with the frequencies Wk given by eq.( [4.16|) 

A detailed analysis and discussion of the choice of initial conditions and the frequen- 
cies ( [4.16|) is provided in the sec. VI. As discussed there, this choice corresponds to the 
large-fc modes being in the conformal adiabatic vacuum state. In what follows we will sub- 
tract the composite operator 'ip'^{t) at the initial time and absorb the term f ('?/'^(^o)) in a 
renormalization of the mass. 

Notice that we have used identical notations in the large and Hartree cases to avoid 
cluttering and also to stress the similarity between the two approximations. In particular, 
we note that the only difference in the expressions for the two cases [eqs.( [4.2^ ) and ( [4.11[ ), 
respectively] is a factor of three appearing in the self interaction term in the equations for 
the zero mode. 

It is instructive to compare the Hartree approximation with the self-consistently one-loop 
approximation presented in sec. II. We see comparing eq.( ^.16D with eq.( [4.27D -( [4.28D that 
they are rather similar. The only difference being that the quantum fluctuations are present 
in both Hartree equations whereas they are only present in the zero mode equation in the 
self-consistently one- loop approximation. In fact, in the weak coupling limit, the Hartree 
approximation becomes the self-consistently one-loop approximation. 



V. RENORMALIZATION 



Renormalization is a very subtle but important issue in gravitational backgrounds [30 
The fluctuation contribution {ip^{x,t)), the energy, and the pressure all need to be renor- 
malized. The renormalization aspects in curved space times have been discussed at length in 



the literature [|30| and have been extended to the large self-consistent approximations for 



the non-equilibrium backreaction problem in P0|, p!^ - p!6| , |4T[] . More recently, a consistent and 
covariant regularization scheme that can be implemented numerically has been proposed 

In terms of the effective mass term for the large limit given by ([4.13| ) and defining the 
quantity 
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Bit)=a\t) [M^{t)-n/Q 



A 



X 



B 



(5.1) 
(5.2) 



where the subscript B stands for bare quantities, we find the following large k behavior for 
the case of an arbitrary scale factor a{t) (with a(0) = 1): 



1 



1 



B{t) 



+ 



d 



3B(ty + a(t)— a(t)B(t) } + 0(l/k 
at ^ -I 



8k^ a2(t) [ 

+ 0{l/k^) , (5.3) 
k 1 
^4(t) + 2ka'^{t) 

+ TTTv^TTTT {-B{t)^ - a{tfB{t) + 3a{t)d{t)B{t) - 4d\t)B{t)} + 0(1/P) 



B{t) + 2a' 



a^{t) 
5« + C(l/fc^) 



(5.4) 



fk{t)m+h{t)n{t) 



d{t) 



ka^{t)a{t) Ak^a'^it) 



2^ Bit) 
a{t) 



+ 0{l/k^ 



(5.5) 



Although the divergences can be dealt with by dimensional regularization, this procedure 
is not well suited to numerical analysis (see however ref . |^ ) . We will make our subtractions 
using an ultraviolet cutoff, Aa(t), constant in physical coordinates. This guarantees that the 
counterterms will be time independent. The renormalization then proceeds much in the 
same manner as in reference |T^; the quadratic divergences renormalize the mass and the 
logarithmic terms renormalize the quartic coupling and the coupling to the Ricci scalar. In 
addition, there is a quartic divergence which renormalizes the cosmological constant while 
the leading renormalizations of Newton's constant and the higher order curvature coupling 
are quadratic and logarithmic respectively. The renormalization conditions on the mass, 
coupling to the Ricci scalar and coupling constant are obtained from the requirement that 
the frequencies that appear in the mode equations are finite [O], i.e: 



m| + ^Bn{t) + ^<p\t) + ^{^\t))s = ml + unit) + ^<p\t) + ^{i^\t))a, (5.6) 

while the renormalizations of Newton's constant, the higher order curvature coupling, and 
the cosmological constant are given by the condition of finiteness of the semi-classical 
Einstein-Friedmann equation: 



(5.7) 



Finally we arrive at the following set of renormalizations 14 1: 
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1 



StiNGr 

Or 

N 
Kr 
N 



1-kNG 



B 



m. 



ln(A/K) 
167r2 



_ / W ln(A//t) 
V^^ 6/ 167r2 



A^ 



IGtt 



2 A2 mi In A//t 
2 ^ 167r2 2 



167r2 



m 



m| + A 



A2 



— Xr C,R — 



ln(A//t) 
167r2 '■"■■"^ 167r2 
1\ ln(A/K) 




M\t) - 



6 



(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 



Here, k is the renormalization point. 

the renormahzations found using dimensional regularization [pT],p 



As expected, the logarithmic terms are consistent with 

Again, we set aji = 



and choose the renormalized cosmological constant such that the vacuum energy is zero in 
the true vacuum. We emphasize that while the regulator we have chosen does not respect 
the covariance of the theory, the renormalized energy momentum tensor defined in this way 
nevertheless retains the property of covariant conservation in the limit when the cutoff is 
taken to infinity. 

The logarithmic subtractions can be neglected because of the coupling A < 10^^^. Using 
the Planck scale as the cutoff and the inflaton mass m^ as a renormalization point, these 
terms are of order Xln[Mpi/m{i] < 10"^'', for m > 10^ GeV . An equivalent statement is 
that for these values of the coupling and inflaton masses, the Landau pole is well beyond the 
physical cutoff Mpi. Our relative error in the numerical analysis is of order 10^^, therefore our 
numerical study is insensitive to the logarithmic corrections. Though these corrections are 
fundamentally important, numerically they can be neglected. Therefore, in the numerical 
computations that follow, we will neglect logarithmic renormalization and subtract only 
quartic and quadratic divergences in the energy and pressure, and quadratic divergences in 
the fluctuation contribution. 



A. Renormalized Equations of Motion for Dynamical Evolution in the Large N limit 

It is convenient to introduce the following dimensionless quantities and deflnitions, 

r = niRt ; h = ; q= ujq = ; g = —- ^ (5.15) 

rriR niR uir Stt^ 

^'M = ^0'W ; 9^ir) = ^{nt))R ; m^V^Mt). (5.16) 

R R 

Choosing C,r = (minimal coupling) and the renormalization point k = \mR\ and setting 
a(0) = 1, the equations of motion become for unbroken symmetry: 
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^ + ^h-^ + 1 + r]\T) + gj:{r) 



dr"^ dr a^ir) 



r]{T) = 



U{r) = 



/,(0) = [-/i(0)-*cu,]/,(0) 



g2 + l + r/2(0) 



6m| 



+ (7S(0) 



We find for broken symmetry, 



r/(r) = , 



-l + r/2 + ^S(r) 



/,(0) = [-/i(O) - z^,] /,(0) 



^-l + r^^(0)-^ + ^S(0) 



g' + l-r/2(0) + 



6m^ 
6m| 



^7S(0) 



for > 1 - r7^(0) + 



for < 1 - 77^(0) + 



6m| 



Here, 



q^dq 



l/.(r)P 



1 



e(g- 1) (m\t) n{T) 



q a[T) 



2g3 



mi 



(5.17) 



(5.18) 



^7S(0) 



-(7S(0). (5.19) 



botli for unbroken and broken symmetry. 

The initial conditions for //(r) will be specified later. An important point to notice is 
that the equation of motion for the g = mode coincides with that of the zero mode ( ^.ISj) . 
Furthermore, for //(r 00) 7^ 0, a stationary (equilibrium) solution of the eq.( p.l(j| ) is 
obtained for broken symmetry when the sum rule |lO,ITl]Jl^,IT3 



- 1+772(00) + ^S(oo) = 



(5.20) 



is fulfilled. This sum rule is nothing but a consequence of Goldstone's theorem and is a 
result of the fact that the large N approximation satisfies the Ward identities associated 
with the 0{N) symmetry, since the term —1 + 77^ + gTi is seen to be the effective mass of 
the modes transverse to the symmetry breaking direction, i.e. the Goldstone modes in the 
broken symmetry phase. 

The renormalized dimensionless evolution equations in the Hartree approximation are 
very similar to eqs.( |5.17| )-( |5.19| ). They can be obtained just dividing by three the 7]'^ term 
in the zero mode equation. [Compare with eqs. ([4.12| )-( [4.13| ) and ( [4.2?| )-( ^^ )]. 
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In terms of the zero mode ?7(r) and the quantum mode function given by eq. (|5.18|) we find 
that the Friedmann equation for the dynamics of the scale factor in dimensionless variables 
is given by 



1 "-0 



2 ^ 4nNm% 



(5.21) 



and the renormalized energy and pressure are given by: 



+ f / <l'dq |/,r - 5«(g, r) + ^- - e(g - 1) 5(^)(g, r 
2Nmi 



(5.22) 



A, 



+ 9 J q'dq [|// - 5«(g, r) + ^ {\fX - e(g - 1) S^'\q, r 



(5.23) 



where the subtractions S^^^ and S^"^^ are given by the right hand sides of eqs. (|5.4]) and ( p3| 
respectively. 

The renormalized energy and pressure are covariantly conserved: 



eR{r) + 3h{T) (p + e)^(r) = 



(5.24) 



From the evolution of the mode functions that determine the quantum fluctuations, we 
can study the growth of correlated domains with the equal time correlation function, 



S{x,t) = {i;{x,t)^{6,t)) 



2(2^ 



Ak-x 



(5.25) 



which can be written in terms of the power spectrum of quantum fluctuations, [/^(t)^. It 
is convenient to define the dimensionless correlation function, 



5(p,r) 



S{\xlt) 



q dq sin[gp] |/g(r) 



p = mR\x\ 



m\ Att'^p Jo 

We now have all the ingredients to study the particular cases of interest. 



(5.26) 



VI. CONFORMAL TIME ANALYSIS AND INITIAL CONDITIONS 

The issue of renormalization and initial conditions is best understood in conformal time 
which is a natural framework for adiabatic renormalization and regularization. 
Quantization in conformal time proceeds by writing the metric element as 



ds^ = C\T){dT^ -dx^). 
Under a conformal rescaling of the field 



(6.1) 
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^x,t)=x{x,T)/C{T), 



(6.2) 



the action for a scalar field (with the obvious generalization to N components) becomes, 
after an integration by parts and dropping a surface term 



with 



S = J d'xdT [^{xT -li^xr -V{X)} 



v(x) = c\r)v{^/c{T)) - c\r)'^x\ 



(6.3) 



(6.4) 



where TZ = 6C"{T)/C^{T) is the Ricci scalar, and primes stand for derivatives with respect 
to conformal time T. 

The conformal time Hamiltonian operator, which is the generator of translations in T, 
is given by 



Ht 



d'x\^Ul+^-iWxr + Vix) 



(6.5) 



with being the canonical momentum conjugate to = x'- Separating the zero mode 

of the field x 



xix,T) = xo{T) + xix,T), 



(6.6) 



and performing the large N or Hartree factorization on the fluctuations we find that the 
Hamiltonian becomes linear plus quadratic in the fluctuations, and similar to a Minkowski 
space-time Hamiltonian with a T dependent mass term given by 



m' + i^-l)n + ^xl{r) + ^{f) 



(6.7) 



We can now follow the steps and use the results of reference [T^ for the conformal time 
evolution of the density matrix by setting a{t) = 1 in the proper equations of that reference 
and replacing the frequencies by 



(6.8) 



and the expectation value in ( |6.7[ ) is obtained in this T evolved density matrix. The time 
evolution of the kernels in the density matrix (see |jl3|) is determined by the mode functions 
that obey 



dT^ 



Fu{T) = 0. 



The Wronskian of these mode functions 

W{F,F*) = FlFt-F,F*^ 



(6.9) 



(6.10) 



is a constant. It is natural to impose initial conditions such that at the initial T the density 
matrix describes a situation of local thermodynamic equilibrium and therefore commutes 
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with the conformal time Hamiltonian at the initial time. This imphes that the initial con- 
ditions of the mode functions Fk(T) be chosen to be (see [ |T3[] ) 



Fk{%) = ; F^(r„) = -t^uJk{%)Fk{%). (6.11) 



With such initial conditions, the Wronskian ( |6.10|) takes the value 

W{F,F*) = -2i. (6.12) 
These initial conditions correspond to the choice of mode functions which coincide with the 



first order adiabatic modes and those of the Bunch-Davies vacuum for large momentum ^0 
To see this clearly, we write the solution of eg. (p.9[) in the form. 



D^(r) = e^T,^'^(^')'^^', (6.13) 

with the function Rk{T) obeying the Riccati equation 

R'f, + Rl + k'' + M''{T) = 0. (6.14) 

This equation posses the solution 

D rn-\ •; I D r'-r\ •-Ri,fc(^) , R2,k{T) .Rz^k{T) . Ri,k{T) , 1 /« i c-n 
Ru{r) = -ik + i?o,fc(^) - I 7 + — p ^ .3 + + ^TF (6-15) 

ixi (\i hj 

and its complex conjugate. We find for the coefficients: 

Ro,k = ; Ri,k = 2-^^(^) ' -^2,A: = —^R^i^k 
1 1 

-R3,fc = 2 {^2,k - Ri,k) ; R4,k = -- (-R's.fc + 2Ri^kR2,k) ■ (6.16) 

The solutions Fk{T) obeying the boundary conditions (|6.11|) are obtained as linear combi- 
nations of this WKB solution and its complex conjugate 

Fk{T) = . \ ^ [(1 + l)D,iT) + (1 - j)DliT)] , (6.17) 



where the coefficient 7 is obtained from the initial conditions. It is straightforward to find 
that the real and imaginary parts are given by 

jR = l + 0{l/k') ■ jj = 0{l/e). (6.18) 



Therefore the \arge-k mode functions satisfy the adiabatic vacuum initial conditions |^ 
This, in fact, is the rationale for the choice of the initial conditions ( 6.11| ). 



Following the analysis presented in we find, in conformal time that 

d^k 
2(27r) 



{x\x,r))= I -^\F,{r)\'. (6.19) 
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The Heisenberg field operators xi^iT) ^ind their canonical momenta n-^(x, T) can be ex- 
panded as: 



n^(f,T) 



V2(27r)3/2 
d^k 
v/2(27r)3/2 



a,F,{T) + al,F;{T) 
a,F',{T) + aUr:' (T) 



Ak-x 



Ak-x 



(6.20) 
(6.21) 



with the time independent creation and annihilation operators at and a], obeying canonical 
commutation relations. Since the fluctuation flelds in comoving and conformal time are 
related by the conformal rescaling (|6.2| ), it is straightforward to see that the mode functions 
in comoving time are related to those in conformal time simply as 



fk{t) 



C{T)- 



(6.22) 



Therefore the initial conditions ( |6.11| ) on the conformal time mode functions imply the initial 
conditions for the mode functions in comoving time are given by eg. ( [4.15[ ). 

For renormalization purposes we need the \axge-k behavior of |/fc(t)P , |/fc(t)p, which are 
determined by the large-/c behavior of the conformal time mode functions and its derivative. 
These are given by 



\F'k(X)\' = k 



1 - 



1 + 



Ri,k{r) 


1 


k^ 




RiA'F) 


1 


k^ 





(6.23) 



We note that the large k behavior of the mode functions to the order needed to renormal- 
ize the quadratic and logarithmic divergences is insensitive to the initial conditions. This is 
not the case when the initial conditions are imposed as described in |T^|T^ . Thus the merit 
in considering the initial conditions in conformal time 



m 



The correspondence with the comoving time mode functions is given by: 

\Fk{r)\' 



\Mt)\' 



1 

CHI) 



\F'k{r)\^ 



H d ' 

c{T)dr^ 



\FkiT)\' 



(6.24) 



There is an important physical consequence of this choice of initial conditions, which is 
revealed by analyzing the evolution of the density matrix. 

In the large N or Hartree (also to one-loop) approximation, the density matrix is Gaus- 
sian, and deflned by a normalization factor, a complex covariance that determines the diag- 
onal matrix elements and a real covariance that determines the mixing in the Schrodinger 
representation as discussed in reference (and references therein). 

In conformal time quantization and in the Schrodinger representation in which the fleld 
X is diagonal the conformal time evolution of the density matrix is via the conformal time 
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Hamiltonian ( |6.5| ). The evolution equations for the covariances is obtained from those given 
in reference [13| by setting a{t) = 1 and using the frequencies ujI{T) = k"^ + Ai'^{T). In 



particular, by setting the covariance of the diagonal elements (given by equation (2.20) in 
IH; see also equation (2.44) of H), 



(6.25) 



we find that with the initial conditions ( |6.11|) , the conformal time density matrix is that 
of local equilibrium at % in the sense that it commutes with the conformal time Hamilto- 
nian. However, it is straightforward to see, that the comoving time density matrix does not 
commute with the comoving time Hamiltonian at the initial time to. 

An important corollary of this analysis and comparison with other initial conditions 
used in comoving time is that assuming initial conditions of local equilibrium in comoving 
time leads to divergences that depend on the initial condition as discussed at length in 
P^ . This dependence of the renormalization counterterms on the initial condition was 
also realized in ref. ||2^ within the context of the CTP formulation. Imposing the initial 
conditions corresponding to local thermal equilibrium in conformal time, we see that: i) 
the renormalization counterterms do not depend on the initial conditions and ii) the mode 
functions are identified with those corresponding to the adiabatic vacuum for large momenta. 
This is why we prefer the initial conditions ( |6.11| ). 

For our main analysis we choose this initial temperature to be zero so that the resulting 
density matrix describes a pure state, which for the large momentum modes coincides with 
the conformal adiabatic vacuum. Such zero temperature choice seems appropriate after the 
exponential inflation of the universe. 

Particle Number: 

We write the Fourier components of the field x ^-^^d its canonical momentum H^ given 
by -{^) as: 



H,,,(T) = i[a,F^(T) + al,F;'(T) 



(6.26) 
(6.27) 



These (conformal time) Heisenberg operators can be written equivalently in terms of the T 
dependent creation and annihilation operators 



Sc^fcC^o) 



nx,fe(T) = -i] 



(6.28) 
(6.29) 



The operators ak{T) ; afc(T) are related by a Bogoliubov transformation. The number 
of particles referred to the initial Fock vacuum of the modes F^, is given by 



N,{T) = {aiiryauiT)) = - 



Fk{T) 
FkiO) 



1 + 



FL{r) 

Fk{T) 



1 

2 ' 



(6.30) 
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or alternatively, in terms of the comoving mode functions fk{t) = Fk{T)/C{T) we find 



Nkit) 



fk{t) 



fkiO) 



fkit) 



Using the large fc-expansion of the conformal mode functions given by eqs. 



(6.31) 



:2l we find 



the large-/c behavior of the particle number to be '^^=^ 0{l/k ), and the total number of 
particles (with reference to the initial state at %) is therefore finite. 



VII. FIELDS EVOLUTION ON A FIXED FRW BACKGROUND 



We consider in this section the matter evolution on radiation or matter dominated FRW 
cosmologies [l^. The case for de Sitter expansion will be discussed in sec. VIII ||T5|| . 

We write the scale factor as a{t) = {t/to)^ with n = 1/2 and n = 2/3 corresponding 
to radiation and matter dominated backgrounds, respectively. Note that the value of to 
determines the initial Hubble constant since 



H{to) 



a (to 



a{to 



n 

To 



We now solve the system of equations ( [4.27] ) - ( [4.29| ) in the Hartree approximation, with 



( [4. lip replacing ( [4.27|) in the large limit. We begin by presenting an early time analysis 
of the slow roll scenario. We then undertake a thorough numerical investigation of various 
cases of interest. For the symmetry broken case, we also provide an investigation of the late 
time behavior of the zero mode and the quantum fluctuations. We use the dimensionless 
variables (|57[5|) -(CT)- 

We will assume minimal coupling to the curvature, C,r = 0. In the cases of interest, 
TZ <^ fi^, so that finite C,r has little effect. 



A. Early Time Solutions for Slow Roll 



For early times in a slow roll scenario [m^ = — yU^, ?7(to) ^ 1]; "we can neglect in eqs. ( [4. 27] ) 
or ( |4.11|) and in eq. ( |4.28| ) both the quadratic and cubic terms in ri{t) as well as the quantum 
fluctuations {ip'^it))^ [recall that ('?/'^(to))r = 0]. Thus, the differential equations for the zero 
mode ( |4.27| ) or ([4.11|) and the mode functions ( |4.28|) become linear equations. In terms 
of the scaled variables introduced above, with a(t) = (n = 2/3 for a matter dominated 
cosmology while n = 1/2 for a radiation dominated cosmology) we have: 
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Uk{t) = . 
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(7.2) 



The solutions to the zero mode equation ( [7.1D are 
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r]{t)=ct-''I,{t)+dt-''K,{t) , 



(7.3) 



where u = {3n — l)/2, and /jy(t) and Ky{t) are modified Bessel functions. Tlie coefficients, c 
and (i, are determined by the initial conditions on rj. For //(to) = ^-^d ?7(to) = 0, we have: 
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(7.5) 



Taking the asymptotic forms of the modified Bessel functions, we find that for intermediate 
times ri{t) grows as 
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(7.6) 



We see that riii) grows very quickly in time, and the approximations ( |7. 1| ) and ( [7. 2] ) will 
quickly break down. For the case shown in fig.l (with n = 2/3, //(to) = 10"'', and ?7(to) = 0), 
we find that this approximation is valid up to t — to — 10. 

The equations for the mode functions ( [7.2| ) can be solved in closed form for the modes 
in the case of a radiation dominated cosmology with n = 1/2. The solutions are 
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Uk{t) = Cfc e-* f/ - - — , -, 2t + 4 e'* M - - — , -,2t] . 
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(7.7) 



Here, U{-) and M(-) are confluent hypergeometric functions ||32| (in another common nota- 
tion, M(-) = iFi(-)), and the and 4 are coefficients determined by the initial conditions 
( [4.30| ) on the modes. The solutions can also be written in terms of parabolic cylinder func- 
tions. 

For large t we have the asymptotic form 



Uk{t) 4 e* (2t)-(3/4+fc^<o/2) 
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(7.8) 



Again, these expressions only apply for intermediate times before the nonlinearities have 
grown significantly. 



B. Numerical Analysis 

We now present the numerical analysis of the dynamical evolution of scalar fields in time 
dependent, matter and radiation dominated cosmo logical backgrounds. We use initial values 
of the Hubble constant such that i/(to) > 0.1. For expansion rates much less than this value 
the evolution will look similar to Minkowski space, which has been studied in great detail 



elsewhere [0,^. As will be seen, the equation of state found numerically is, in the majority 
of cases, that of cold matter. We therefore use matter dominated expansion for the evolution 
in much of the analysis that follows. While it presents some inconsistency at late times, the 
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evolution in radiation dominated universes remains largely unchanged, although there is 
greater initial growth of quantum fluctuations due to the scale factor growing more slowly 
in time. Using the large and Hartree approximations to study theories with continuous 
and discrete symmetries respectively, we treat three important cases. They are 1) < 0, 
r]{to) < 1; 2) < 0, r/(to) > 1; 3) > 0, r]{to) > 1. 

In presenting the figures, we have shifted the origin of time such that t t' = t — to. 
This places the initial time, to, at the origin. In these shifted coordinates, the scale factor 
is given by 



where, once again, n = 2/3 and n = 1/2 in matter and radiation dominated backgrounds 
respectively, and the value of r is determined by the Hubble constant at the initial time: 



Case 1: < 0, i]{to) ^ 1. This is the case of an early universe phase transition in 
which there is little or no biasing in the initial configuration (by biasing we mean that the 
initial conditions break the t] — >■ —t] symmetry). The transition occurs from an initial tem- 
perature above the critical temperature, T > Tc, which is quenched at to to the temperature 
Tf <^ Tc. This change in temperature due to the rapid expansion of the universe is modeled 
here by an instantaneous change in the mass from an initial value = T^/T^^ — 1 to a final 
value mj = — 1. We will use the value = 1 in what follows. This quench approximation 
is necessary since the low momentum frequencies ([4.16|) appearing in our initial conditions 
( [4.3U| ) are complex for negative mass squared and small ri{to). An alternative choice is to 
use initial frequencies given by 



These frequencies have the attractive feature that they match the conformal adiabatic fre- 
quencies given by ( |4.16|) for large values of k while remaining positive for small k. We find 
that such a choice of initial conditions changes the quantitative value of the particle number 
by a few percent, but leaves the qualitative results unchanged. 

We plot the the zero mode fjit), the equal time correlator gTj{t), the total number of 
produced particles gN(t) (see sec. VI for a discussion of our definition of particles), the 
number of particles gNk{t) as a function of wavenumber for both intermediate and late 
times, and the ratio of the pressure and energy densities p(t)/e(t) (giving the equation of 
state). 

Figs, la-e shows these quantities in the large approximation for a matter dominated 
cosmology with an initial condition on the zero mode given by ri{to = 0) = 10~^, ?7(to = 0) = 
and for an initial expansion rate of H{to) = 0.1. This choice for the initial value of rj stems 
from the fact that the quantum fluctuations only have time to grow significantly for initial 
values satisfying rjito) <C y/g; for values rjito) ^ y/g the evolution is essentially classical. 
This result is clear from the intermediate time dependence of the zero mode and the low 
momentum mode functions given by the expressions (|7.6|) and ( [7.8|) respectively. 
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FIGURES 
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FIG. 1. Symmetry broken, slow roll, large N, matter dominated evolution of (a) the zero mode 
r/(t) vs. t, (b) the quantum fluctuation operator g'S{t) vs. t, (c) the number of particles gN{t) 
vs. t, (d) the particle distribution gNk{t) vs. k &i t = 149.1 (dashed line) and t = 398.2 (solid 
line), and (e) the ratio of the pressure and energy density p{t)/e(t) vs. t for the parameter values 
m2 = -1, r7(to) = 10-^ f]{to) = Q,g = IQ-^^, H{to) = 0.1. 

After the initial growth of the fluctuation gT,{t) (fig. lb) we see that the zero mode (fig. la) 
approaches the value given by the minimum of the tree level potential, = 1, while gT^{t) 
decays for late times as 

^^('^ " ^ = ^ • 

For these late times, the Ward identity corresponding to the 0{N) symmetry of the field 
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theory is satisfied, enforcing tlie condition 



-l + n\t)+gnt) = Q. (7.9) 
Hence, tfie zero mode approacfies tfie classical minimum as 

Figure Ic depicts the number of particles produced. After an initial burst of particle 
production, the number of particles settles down to a relatively constant value. Notice that 
the number of particles produced is approximately of order 1/g. In fig. Id, we show the 
number of particles as a function of the wavenumber, k. For intermediate times we see the 
simple structure depicted by the dashed line in the figure, while for late times this quantity 
becomes concentrated more at low values of the momentum k. 

Finally, fig.le shows that the field begins with a de Sitter equation of state p = —e but 
evolves quickly to a state dominated by ordinary matter, with an equation of state (averaged 
over the oscillation timescale) p = 0. This last result is a bit surprising as one expects from 
the condition ( [7.9|) that the particles produced in the final state are massless Goldstone 
bosons which should have the equation of state of radiation. However, as shown in fig. Id, 
the produced particles are of low momentum, g -C 1, and while the effective mass of the 
particles is zero to very high accuracy when averaged over the oscillation timescale, the 
effective mass makes small oscillations about zero so that the dispersion relation for these 
particles differs from that of radiation. In addition, since the produced particles have little 
energy, the contribution to the energy density from the zero mode, which contributes to a 
cold matter equation of state, remains significant. 

In figs. 2a-e we show the same situation depicted in fig.l using the Hartree approximation. 
The initial condition on the zero mode is rjitQ = 0) = -\/3 ■ 10"''; the factor of -\/3 appears 
due to the different scaling in the zero mode equations, ( |4.27] ) and ( |5.18| ), which causes the 



minimum of the tree level effective potential in the Hartree approximation to have a value 
oirf = v^. Again, the Hubble constant has the value H{to) =0.1. Here, we see again that 
there is an initial burst of particle production as gTi{t) (fig. 2b) grows large. However, the 
zero mode (fig. 2a) quickly reaches the minimum of the potential and the condition 

-l+r]\t)/3 + g^{t) = (7.10) 

is approximately satisfied by forcing the value of gT,(t) quickly to zero. 

There are somewhat fewer particles produced here compared to the large N case, and 
the distribution of particles is more extended. Since the effective mass of the particles 
is nonzero, we expect a matter dominated equation of state (fig 2e) for later times. The 
fact that the Hartree approximation does not satisfy Goldstone's theorem means that the 
resulting particles must be massive, explaining why somewhat fewer particles are produced. 

Finally, we show the special case in which there is no initial biasing in the field, 77(^0 = 
0) = 0, fi{tQ = 0) = 0, and H{to) = 0.1 in figs. 3a-d. With such an initial condition, the 
Hartree approximation and the large limit are equivalent. The zero mode remains zero 
for all time, so that the quantity gT^it) (fig. 3a) satisfies the sum rule (|7.9|) by reaching the 
value one without decaying for late times. Notice that many more particles are produced in 
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this case (fig 3b); the growth of the particle number for late times is due to the expansion 
of the imiverse. The particle distribution (fig.3c) is similar to that of the slow roll case in 
fig.l. The equation of state (fig.3d) is likewise similar. 
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FIG. 2. Symmetry broken, slow roll, Hartree, matter dominated evolution of (a) the zero mode 
rj{t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles gN(t) 
vs. t, (d) the particle distribution gNk{t) vs. k at t = 150.7 (dashed line) and t = 396.1 (solid 
line), and (e) the ratio of the pressure and energy density p{t)le{t) vs. t for the parameter values 
m2 = -1, r?(to) = . 10-7, r]{to) = 0, 5 = IQ-^^, H{to) = 0.1. 

In each of these cases of slow roll dynamics, increasing the Hubble constant has the 
effect of slowing the growth of both 77 and g'^{t). The equation of state will be that of a de 
Sitter universe for a longer period before moving to a matter dominated equation of state. 
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Otherwise, the dynamics is much the same as in figs. 1-3. 

Case 2: < 0, r]{tQ) S> 1. We now examine the case of a chaotic infiationary scenario 
with a symmetry broken potential. In chaotic infiation, the zero mode begins with a value 
rj{t) ^ 1. During the de Sitter phase, H ^ 1, and the field initially evolves classically, 
dominated by the first order derivative term appearing in the zero mode equation (see 
( [4.271 ) and ( [4.28|) ). Eventually, the zero mode rolls down the potential, ending the de Sitter 
phase and beginning the FRW phase. We consider the field dynamics in the FRW universe 
after the end of infiation. We thus take the initial temperature to be zero, T = 0. 
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FIG. 3. Symmetry broken, no roll, matter dominated evolution of (a) the quantum fluctuation 
operator gT,{t) vs. t, (b) the number of particles gN{t) vs. t, (c) the particle distribution gNk{t) 
vs. k at t = 150.1 (dashed line) and t = 397.1 (solid line), and (d) the ratio of the pressure and 
energy density p{t)/e{t) vs. t for the parameter values = —1, ry(to) = 0, ?7(to) = 0, = 10~^^, 
H{to) = 0.1. 

Figure 4 shows our results for the quantities, ri{t), gT,{t), gN{t), gNk{t), and p{t)/e{t) 
for the evolution in the large approximation within a radiation dominated gravitational 
background with H{tQ) = 0.1. The initial condition on the zero mode is chosen to have the 
representative value ri{to = 0) = 4 with fi{tQ = 0) = 0. Initial values of the zero mode much 
smaller than this will not produce significant growth of quantum fiuctuations; initial values 
larger than this produces qualitatively similar results, although the resulting number of 
particles will be greater and the time it takes for the zero mode to settle into its asymptotic 
state will be longer. 
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FIG. 4. Symmetry broken, chaotic, large A^, radiation dominated evolution of (a) the zero mode 
r/(t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles gN{t) vs. 
t, (d) the particle distribution gNk{t) vs. k at t = 76.4 (dashed line) and t = 392.8 (solid line), and 
(e) the ratio of the pressure and energy density p{t)/e{t) vs. t for the parameter values = — 1, 
r?(to) = 4, i){to) = 0,g = W'^^ H{to) = 0.1. 

We see from fig. 4a that the zero mode oscillates rapidly, while the amplitude of the os- 
cillation decreases due to the expansion of the universe. This oscillation induces particle 
production through the process of parametric amplification (fig.4c) and causes the fluc- 
tuation 5'S(t) to grow (fig.4b). Eventually, the zero mode loses enough energy that it is 
restricted to one of the two minima of the tree level effective potential. The subsequent 
evolution closely follows that of Case 1 above with gTj{t) decaying in time as l/a^(t) ~ 1/t 
with T] given by the sum rule (|7.9| ). The spectrum (fig.4d) indicates a single unstable band 
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of particle production dominated by the modes /c = 1/2 to about k = 3 for late times. The 
structure within this band becomes more complex with time and shifts somewhat toward 



lower momentum modes. Such a shift is also observed in Minkowski spacetimes |T0|-p!2 



Figure 4e shows the equation of state which we see to be somewhere between the relations 
for matter and radiation for times out as far as t = 400, but slowly moving to a matter 
equation of state. Since matter redshifts as l/a^{t) while radiation redshifts as l/a^(t), the 
equation of state should eventually become matter dominated. Given the equation of state 
indicated by fig.4e, we estimate that this occurs for times of order t = 10^. The reason 
the equation of state in this case differs from that of cold matter as was seen in figs. 1-3 is 
that the particle distribution produced by parametric amplification is concentrated at higher 
momenta, A; ~ 1. 

Figure 5 shows the corresponding case with a matter dominated background. The results 
are qualitatively very similar to those described for fig. 4 above. Due to the faster expansion, 
the zero mode (fig. 5a) finds one of the two wells more quickly and slightly less particles are 
produced. For late times, the fluctuation gT,{t) (fig. 5b) decays as l/a^{t) oc l/f^^^. Again we 
see an equation of state (figs. 5e) which evolves from a state between that of pure radiation 
or matter toward one of cold matter. 

The Hartree case is depicted in fig. 6 for a matter dominated universe, with the initial 
condition on the zero mode 7^(^0 = 0) = 4\/3. Again, the evolution begins in much the same 
manner as in the large approximation with oscillation of the zero mode (fig. 6a), which 
eventually settles into one of the two minima of the effective potential. Whereas in the 
large A^ approximation, the zero mode approaches the minimum asymptotically [as given by 
( [7.9| ) and our late time analysis below], in the Hartree approximation we see that the zero 
mode finds the minimum quickly and proceeds to oscillate about that value. The two point 
correlator (fig.6b) quickly evolves toward zero without growing large. Particle production 
in the Hartree approximation (figs. 6c-d) is again seen to be inefficient compared to that of 
the large A^ case above. Fig. 6e again shows that the equation of state is matter dominated 
for all but the earliest times. 

A larger Hubble constant prevents significant particle production unless the initial ampli- 
tude of the zero mode is likewise increased such that the relation rjito) ^ H(to) is satisfied. 
For very large amplitude T]{to) ^ 1, to the extent that the mass term can be neglected 
and while the quantum fluctuation term has not grown to be large, the equations of motion 
( [4.271 ), (|4.28|) , and (|4.11|) are scale invariant with the scaling t] — >■ ^t], H fiH, t t/fi, 
and k ^ ^k, where /i is an arbitrary scale. 

Case 3: > 0, rjito) ^ 1. The final case we examine is that of a simple chaotic 
scenario with a positive mass term in the Lagrangian. Again, the FRW stage occurs after 
the inflationary expansion; this allows us to take zero initial temperature. 

Figure 7 shows this situation in the large A^ approximation for a matter dominated 
cosmology. The zero mode, ri{t), oscillates in time while decaying in amplitude from its 
initial value of 77(^0 = 0) = 5, r](to = 0) = (fig. 7a), while the quantum fluctuation, 5'E(i), 
grows rapidly for early times due to parametric resonance (figs. 7b). We choose here an 
initial condition on the zero mode which differs from that of figs 4-5 above since there is no 
significant growth of quantum fluctuations for smaller initial values. From fig.7d, we see 
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FIG. 5. Symmetry broken, chaotic, large N, matter dominated evolution of (a) the zero mode 
ri{t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles gN{t) vs. 
t, (d) the particle distribution gNk{t) vs. k aX t = 50.8 (dashed line) and t = 399.4 (solid line), and 
(e) the ratio of the pressure and energy density p{t)/e{t) vs. t for the parameter values = — 1, 
r?(to) = 4, riito) = 0, g = IQ-'^ H{to) = 0.1. 

that there exists a single unstable band at values of roughly A; = 1 to A; = 3, although 
careful examination reveals that the unstable band extends all the way to A; = 0. The 
equation of state is depicted by the quantity p{t)/e{t) in fig.7e. As expected in this massive 
theory, the equation of state is matter dominated. 



34 




FIG. 6. Symmetry broken, chaotic, Hartree, matter dominated evolution of (a) the zero mode 
rj{t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles gN(t) 
vs. t, (d) the particle distribution gNk{t) vs. k at t = 151.3 (dashed line) and t = 397.0 (solid 
line), and (e) the ratio of the pressure and energy density p{t)le{t) vs. t for the parameter values 
m? = -1, rj{to) = 4 ■ 3V2, ri{to) = 0,g = IQ-^^, H{to) = 0.1. 

The final case is the Hartree approximation, shown in fig. 8. Here, parametric amplifi- 
cation is entirely inefficient when expansion of the universe is included and we require an 
initial condition on the zero mode of ri{to = 0) = 12^/3 to provide even meager growth of 
quantum fluctuations. We have used a matter dominated gravitational background with 
if (to) =0.1. We see that while the zero mode oscillates (flg.Sa), there is little growth in 
quantum fluctuations (flg.Sb) and few particles produced (flg.Sc). Examining the particle 
distribution (flg.Sd), it is found that the bulk of these particles is produced within a single 
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FIG. 7. Symmetry unbroken, chaotic, large A^, matter dominated evolution of (a) the zero 
mode ri{t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles 
gN{t) vs. t, (d) the particle distribution gNk{t) vs. k att = 17 A (dashed line) and t = 399.7 (solid 
line), and (e) the ratio of the pressure and energy density p{t)le{t) vs. t for the parameter values 



m 



+1, r?(to) = 5, 77(to) = 0, 5 = 10-12, H{tQ) = 0.1. 



resonance band extending from A;~15toA;~16. This resonance develops at early time 
during the large amplitude oscillation of the zero mode. These results are explained by a 
simple resonance band analysis described below. 

At first glance, it is not entirely clear why there are so many more particles produced 
in the large N case of fig. 7 than in the Hartree case of fig. 8. Since in the present case the 
Hubble time is long compared to the oscillation timescale of the zero mode, if <C 1, we 
would expect a forbidden band for early times at the location given approximately by the 
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Figure 8a 
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Figure 8b 
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100 200 300 400 



10 15 20 



Figure 8e 
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FIG. 8. Symmetry unbroken, chaotic, Hartree, matter dominated evolution of (a) the zero 
mode r/(t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles 
gN{t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 50.5 (dashed line) and t = 391.2 (solid 
line), and (e) the ratio of the pressure and energy density p[t)/e{t) vs. t for the parameter values 



m 



+1, r?(to) = 12 • 3I/2, ^(to) = 0, ^7 = 10-12, F(to) = 0.1. 



Minkowski results provided in Ref. [0. In fact, we find this to be the case. 

The solution to this problem is inherent in the band structure of the two cases when 
combined with an understanding of the dynamics in an expanding spacetime. First, we 
note that, for early times when gTj{t) -C 1, the zero mode is well fit by the function ri{t) = 
Vof{t)/O'{t) where f{t) is an oscillatory function taking on values from —1 to 1. This is 
clearly seen from the envelope function ?7o/a(t) shown in fig. 8a (recall that 5'S(t) -C 1 during 
the entire evolution in this case). Second, the momentum that appears in the equations 
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for the modes ( |4.28|) is the physical momentum k/a{t). We therefore write the approximate 



expressions for the locations of the forbidden bands in FRW by using the Minkowski results of 
ri| with the substitutions tiq '-jtiq/ a^{t) (where the factor of 7 accounts for the difference 



in the definition of the non-linear coupling between this study and [|rT|) and — k"^ /a?{t). 

Making these substitutions, we find for the location in comoving momentum k of the 
forbidden band in the large N (fig. 7) and Hartree (fig. 8) cases: 

2 

< k^< y, (large A^) (7.11) 



^ + 3a2(t)<A;2< a\t) 



''^ +47i+4 + ll , (Hartree). (7.12) 



2 ' ' - - " ^ ^ 3a4(t) a2(t) 

The important feature to notice is that while the location of the unstable band (to a first 
approximation) in the case of the continuous 0{N) theory is the same as in Minkowski and 
does not change in time, the location of the band is time dependent in the discrete theory 
described by the non-perturbative Hartree approximation. 
While r]o/a(t) ^ 1, the Hartree relation reduces to 

I < P < ^. (7.13) 

This is the same as the Minkowski result for large amplitude, and one finds that this expres- 
sion accurately predicts the location of the resonance band of fig.8d. However, with time 
the band shifts its location toward higher values of comoving momentum as given by ( [7. 121) , 
cutting off particle production in that initial band. There is continuing particle production 
for higher modes, but since the Floquet index is decreased due to the reduced amplitude 
of the zero mode, since there is no enhancement of production of particles in these modes 
(as these modes begin with at most of order 1 particles), and because the band continues 
to shift to higher momenta while becoming smaller in width, this particle production never 
becomes significant. 

As in the symmetry broken case of figs. 4-6, the equations of motion for large amplitude 
and relatively early times are approximately scale invariant. In fig. 9 we show the case of 
the large N evolution in a radiation dominated universe with initial Hubble constant of 
H[to) = 2 with appropriately scaled initial value of the zero mode of ^^(to) = 16. Again, the 
qualitative dynamics remains largely unchanged from the case of a smaller Hubble constant. 



C. Late Time Behavior 



We see clearly from the numerical evolution that in the case of a symmetry broken 
potential, the late time large solutions obey the sum rule ( |5.20| ). This sum rule is a 
consequence of the late time Ward identities which enforce Goldstone's Theorem. Because 
of this sum rule, we can write down the analytical expressions for the late time behavior of 
the fluctuations and the zero mode. Using (|5.20|) , the mode equation ( [4.28|) becomes 



d{t) d P 



Uk{t) = 0. 



(7.14) 
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This equation can be solved exactly if we assume a power law dependence for the scale factor 
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FIG. 9. Symmetry unbroken, chaotic, large N, radiation dominated evolution of (a) the zero 

mode r]{t) vs. t, (b) the quantum fluctuation operator gT,{t) vs. t, (c) the number of particles 
gN{t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 102.1 (dashed line) and t = 251.6 
(solid line), and (e) the ratio of the pressure and energy density p{t)/£{t) vs. t for the parameter 
values = +1, ??(to) = 16, ?7(to) = 0,9 = 10-^^ H{to) = 2.0. 

a{t) = (t/io)"- The solution is 

(l.fnfl-n\ / h.4.n4.l-n\ 

n-l ) 2-2" \n — l) 

where Jy and are Bessel and Neumann functions respectively, and the constants and 
dfc carry dependence on the initial conditions and on the dynamics up to the point at which 
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the sum rule is satisfied. 

These functions have several important properties. In particular, in radiation or matter 
dominated universes, n < 1, and for values of wavenumber satisfying k 3> t"*-^""-*/^?)) 
mode functions decay in time as l/a{t) ~ Since the sum rule applies for late times, 
t — to ^ 1 in dimensionless units, we see that all values of k except a very small band about 
k = redshift as l/a{t). The k = mode, however, remains constant in time, explaining 
the support evidenced in the numerical results for values of small k (see figs. 1,3). These 
results mean that the quantum fluctuation has a late time dependence of {ip'^{t))r ~ l/a^{t). 
The late time dependence of the zero mode is given by this expression combined with the 
sum rule ( [5.20[ ). These results are accurately reproduced by our numerical analysis. Note 
that qualitatively this late time dependence is independent of the choice of initial conditions 
for the zero mode, except that there is no growth of modes near k = in the case in which 
particles are produced via parametric amplification (figs. 4,5). 

For the radiation n = 
elementary functions: 



and matter dominated n = universes, eq.( |7.15|) reduces to 



a{t) Uuit) = Cfc e^<'"^"' + 4 e-^<'"''^" (RD) 



a(t) Uu{t) = cu e'<'''''' 



1 + 



3ktTt^/-^ 



+ dke ^ 



1 - 



(MD) . (7.16) 



It is also of interest to examine the n > 1 case. Here, the modes of interest satisfy the 
condition k <C t"'~^/tQ for late times. These modes are constant in time and one sees that 
the modes are frozen. In the case of a de Sitter universe, we can formally take the limit 
n oo and we see that all modes become frozen at late times. This case is detailed in sec. 

VIII ra. 



D. Discussion and Conclusions for the FRW background 

We have shown that there can be significant particle production through quantum fluctu- 



ations after inflation [|1^. However, this production is somewhat sensitive to the expansion 
of the universe. From our analysis of the equation of state, we see that the late time dy- 
namics is given by a matter dominated cosmology. We have also shown that the quantum 
fluctuations of the inflaton decay for late times as l/a^(t), while in the case of a symmetry 
broken inflationary model, the inflaton field moves to the minimum of its tree level poten- 
tial. The exception to this behavior is the case when the inflaton begins exactly at the 
unstable extremum of its potential for which the fluctuations grow out to the minimum of 
the potential and do not decay. Initial production of particles due to parametric amplifi- 
cation is significantly greater in chaotic scenarios with symmetry broken potentials than in 
the corresponding theories with positive mass terms in the Lagrangian, given similar initial 
conditions on the zero mode of the inflaton. 
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VIII. FIELDS EVOLUTION ON A FIXED INFLATIONARY BACKGROUND 

(THE DE SITTER UNIVERSE) 



We describe in this section the matter evolution on a fixed de Sitter background using 
large N and Hartree approaches. The evolution with a dynamical background is treated in 
sec. IX. 

A. Evolution for ^(0) = ^(0) = 0. Analytical Results 

We begin by considering the broken symmetry situation in which the expectation value 
of the inflaton field sits atop the potential hill with zero initial velocity. This situation 
is expected to arise if the system is initially in local thermodynamic equilibrium an ini- 
tial temperature larger than the critical temperature and cools down through the critical 
temperature in the absence of an external field or bias. 

The order parameter and its time derivative vanish in the local equilibrium high tem- 
perature phase, and this condition is a fixed point of the evolution equation for the zero 
mode of the infiaton. There is no rolling of the infiaton zero mode in this case, although the 
fiuctuations will grow and will be responsible for the dynamics. 

We can understand the early stages of the dynamics analytically as follows. For very 
weak coupling and early time we can neglect the backreation in the mode equations, which 
in both the large and Hartree cases become. 



— + 3/.- + ^ 



- 1 



(8.1) 



/q(0) = —— ; fgir) = -iy/uJ^ ; ujg = ^Jq^ + r"^-! . 



(8.2) 



The solutions are of the form. 



3 q / 1 9 

fq{T) = exp[--/iT] {a{q) J^{z) + h{q) J-^{z)} ; z ^ - exp[-/iT] ; i/ = W — + -, (8.3) 

where the coefficients a{q) and b{q) are determined by the initial conditions: 



2h sin i/TT 



h 



" \h 



(8.4) 



a{q) 



TT q 



2h sin UTT 



(8.5) 



For long times, e > q/h, these mode functions grow exponentially. 



/,(r) ^ b{q) J.,{z) 



b{q) 



r(i -u) \q 



(8.6) 
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The Bessel functions appearing in the expression for the modes fqir) can be approxi- 
mated by their series expansion, 



LUr, 



4/^2 



h 



^(u-?,/2)hT 



17) 



This is an expansion in powers of /{vh"^) and we conclude that 5'S(r) is dominated by the 
modes with q < \/h. 

The integral for gTi(T) can be approximated by keeping only the modes q < fVh, 
where / is a number of order one, and by neglecting the subtraction term which will cancel 
the contributions from high momenta. Numerically, even with the backreaction taken into 
account, the integral is dominated by modes g < / ^ 10 — 20 in all of the cases that we 
studied (see fig. 12 below). 

The contribution to the fluctuations from these unstable modes is: 



rm 



R 



27r 



1 + 



M2 



^{2u-3)hT 



rrir 



(8.8) 



where again, we have taken the high temperature limit, Ti niR- 

From this equation, we can estimate the value of r^, the 'spinodal time', at which the 
contribution of the quantum fluctuations becomes comparable to the contribution from the 
tree level terms in the equations of motion. This time scale is obtained from the condition 
g^rs) = 0(1): 



1 



{2u-3)h 



In 



3 27: ttirM^ Z 



m 



R/ J 



(8.9) 



which is in good agreement with our numerical results, as will become clear below (see figs. 
10, 15 and 18). For values of h > 1, which, as argued below, lead to the most interesting 
case, an estimate for the spinodal time is. 



3h 



ln[l/V^ + 0(l) 



(8.10) 



which is consistent with our numerical results (see fig. 10). 

For r > Ts, the effects of backreaction become very important, and the contribution from 
the quantum fluctuations competes with the tree level terms in the equations of motion, 
shutting-off the instabilities. Beyond r^, only a full numerical analysis will capture the 
correct dynamics. 

It is worth mentioning that had we chosen zero temperature initial conditions, then the 
coupling g ^ g (see ([T9|)) and the estimate for the spinodal time would have been. 



3h 

r,^-\n[l/g] + 0{l), 



(8.11) 



that is, roughly a factor 2 larger than the estimate for which the de Sitter stage began at 
a temperature above the critical value. Therefore (|8.10|) represents an underestimate of the 
spinodal time scale at which fluctuations become comparable to tree level contributions. 
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The number of e-folds occurring during the stage of growth of spinodal fluctuations is 
therefore, 



or in the zero temperature case, 



Me^—Hl/V^l (8.12) 



M.^''-^Hl/gl (8.13) 



which is a factor 2 larger. Thus, it becomes clear that with g ^ 10^^^ and h >2, a required 
number of e-folds, A/'e ~ 100 can easily be accommodated before the fluctuations become 
large, modifying the dynamics and the equation of state. 

The implications of these estimates are important. The first conclusion drawn from these 
estimates is that a 'quench' approximation is well justified (see fig. 10). While the temper- 
ature drops from an initial value of a few times the critical temperature to below critical 
in just a few e-folds, the contribution of the quantum fluctuations needs a large number of 
e-folds to grow to compensate for the tree-level terms and overcome the instabilities. Only 
for a strongly coupled theory is the time scale for the quantum fluctuations to grow short 
enough to restore local thermodynamic equilibrium during the transition. 

The second conclusion is that most of the growth of spinodal fluctuations occurs during 
the inflationary stage, and with g ^ 10~^^ and H > ttir, the quantum fluctuations become 
of the order of the tree-level contributions to the equations of motion within the number of 
e-folds necessary to solve the horizon and flatness problems. Since the fluctuations grow to 
become of the order of the tree level contributions at times of the order of this time scale, 
for larger times they will modify the equation of state substantially and will be shown to 
provide a graceful exit from the inflationary phase within an acceptable number of e-folds. 

For T < Ts, when the contribution from the renormalized quantum fluctuations can be 
ignored, the Hubble constant is given by the classical contribution to the energy density. In 



terms of the dimensionless quantities introduced above ( |5.16| ), we have 



2 4^' ; 



(8.14) 



In the situation we consider here, with i) = t] = 0, the condition that h > 2 for g ~ 10~^^ 
translates into ~ 10^^ GeV , which is an acceptable bound on the inflaton mass. 

To understand more clearly whether or not the effect of quantum fluctuations and growth 
of unstable modes during the inflationary phase transition can provide a graceful exit sce- 
nario, we must study in detail the contribution to the energy and the equation of state of 
these quantum fluctuations. 

Although we are working in a fixed de Sitter background, the energy and pressure will 
evolve dynamically. A measure of the backreaction effects of quantum fluctuations on the 
dynamics of the scale factor is obtained from defining the 'effective Hubble constant', 

n^ir)^^Hr). (8.15) 
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Therefore, the quantities, 



n{T) 



n{0) \ e{0) 



sir) 



and 



1 + 



p(r) 



give dynamical information of the effects of the backreaction of the quantum fluctuations 
on the dynamics of the scale factor. Whenever p(r) + e^r) ^ 0, Tiir) /TiiQi) ^ 1, or 
'H{t) /1H?{t) 7^ 0, the backreaction from the quantum fluctuations will dramatically change 
the dynamics of the scale factor, and it will no longer be consistent to treat the scale factor 
as fixed. When TC{t) /TC{0) -C 1 the de Sitter era will end. 

From this point onwards only a full treatment of the backreaction, including the correct 
dynamics of the scale factor, will describe the physics. This will be the subject of sec. IX 

m 



B. Numerical Analysis 



We now solve numerically the large N set of equations ( p. 181 ) with the initial conditions 
( Fl9| ), taking r]{0) = r/(0) = 0. 

The numerical code is based on a fourth order Runge-Kutta algorithm for the differential 
equation and an 11-points Newton-Cotes algorithm for the integral, with a typical relative 
errors 10~^ in the differential equation and in the integrals. We have tested for cutoff 
insensitivity with cutoffs Qmax ~ 50 ; 100 ; 150 with no appreciable variation in the numerical 
results. The reason for this cutoff insensitivity is due to the fact that only long- wavelength 
modes grow in amplitude to become non-perturbatively large, whereas the short-wavelength 
modes always have perturbatively small amplitudes. We have chosen r = Tj/T^ = 2 as a 
representative value and g = 10~^^. 

As argued previously, for g ^ 10^^^, the cosmologically interesting time scales for the 
spinodal instabilities to grow during say the first 60-100 e-folds of inflation occur for h> 1, 
leading to H > itlr > 10^'^ GeV which is a phenomenologically acceptable range for the 
Hubble constant during the inflationary stage. 

Fig. 10 shows the contribution from the quantum fluctuations, 5'S(r) vs. r for g = 
10^^^ ; Tj/Tc = 2 ; h = 2 ; //(O) = ; ?7(0) = 0. The quantum fluctuations, as measured 
by (7S(r), grow to be of order 1 in a time scale r ~ 40 which is the time scale predicted by 
the early time estimate ( ^.10| ). We have checked numerically that the energy is covariantly 



conserved, obeying the relation e+3H{p+e) = to our numerical accuracy of one part in 10''. 
Fig. 11 shows H{t) /H{0) vs. r. This figure shows clearly that when the spinodal quantum 
fluctuations become comparable to the tree level contribution to the equations of motion, 
the backreaction on the scale factor becomes fairly large. At this point, the approximation 
of keeping a fixed background breaks down and the full self-consistent dynamics will have to 
be studied. At this time, the inflationary stage basically ends since Ti. is no longer constant. 
This occurs for t ^ 40 giving about 80 e-folds of inflation during the time in which Ti 
is approximately constant and equal to H. Therefore, this new mechanism of spinodal 
fluctuations, with the zero mode sitting atop the potential hill provides a graceful exit of 
the inflationary era without any further assumptions on the evolution of the scalar field. 
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FIG. 10. c/S(r) vs. T for r/(0) = 77(0) = 0; A = lO'^^. ^^2; h = 2. 

Figure Ic 




FIG. 11. H{t)/H{0) vs t for the same parameters as in fig.lO. 



Figure le 




FIG. 12. gS{q,T) vs. g for r = 60. Same parameters as in fig.lO. 



These fluctuations translate into an amplification of the power spectrum at long wave- 
lengths for q ^ h. To see this clearly we plot gS{q,T) = 5'|/g(T)p vs. g for r = 60 in 
fig. 12. This quantity is very small, because of the coupling constant in front, but for r ~ 
it grows to be of order one for long wavelengths (see also fig. 14) and vanishes very fast for 
g > 10. The integral in g'S(r) is dominated by these long wavelengths that become non- 
perturbatively large, whereas the contribution from the short wavelengths remains always 
perturbativcly small. This is the justification for the approximations performed early that 
involved only the long-wavelength modes and cutoffs of order Vh. The equal time spatial 
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correlation function given by eq.( |5.26|) can now be computed explicitly. Figure 13 shows 
S{p;t)/S{0;t) as a function of p for r > 2. We define the correlation length ^(r) as the 
value of p for which the ratio is 1/e. Figure (l.g) shows ^(r); notice that the correlation 
length saturates to a value ^(oo) ~ and that the correlated regions are of horizon size. 

We have performed numerical analysis varying h with the same values of g and for the 
same initial conditions, and found that the only quantitative change is in the time scale for 
(y'E(r) to be of order one. We find that the spinodal time scale grows almost linearly with h 
and its numerical value is accurately described by the estimate ( |8.10| ). The case in which the 
Hubble constant is /i = 0.1 is shown explicitly. Figure 15 shows gJ^lr), which demonstrates 
the oscillatory behavior similar to what is seen in Minkowski space |iy,[TT|. The correlation 
length, ^(t), is shown in fig. 16; its asymptotic value is again approximately given by 1/h. 




FIG. 13. S{p,t)/S{0,t) vs. /? for t > 2. Same parameters as in fig. 10. 

Figure Ih 




FIG. 14. In [|/g(T)|^] vs. r for q = 0,A, 10. Same parameters as in fig. 10. 
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FIG. 15. 5E(t) vs. r for ; r/(0) = 77(0) = 0; ; A = lO'^^ r = 2; h = 0.1. 

Figure 2b 
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FIG. 16. .^(r) vs r. Same parameters as in fig. 15. 



C. The late time limit 



For times r > ^ 40 (for the values of the parameters used in figs. (1)) we see from figs. 
10-14 that the dynamics freezes out. The fluctuation, 5'S(r) = 1, and the mode functions 
effectively describe free, minimally coupled, massless particles. The sum rule. 



1 + ^S(oo) = 0, 



(8.16) 



is obeyed exactly in the large N limit as in the Minkowski case [jl0|-[T^ . 

For the Hartree case g — *■ 3g, but the physical phenomena are the same, with the only 
difference that the sum rule now becomes gTi{oo) = 1/3. We now show that this value is 
a self-consistent solution of the equations of motion for the mode functions, and the only 
stationary solution for asymptotically long times. 

In the late time limit, the effective time dependent mass term, — 1 + 77^ + gT,, in the 
equation for the mode functions, ( 5.18|) , vanishes (in this case with 77 = 0). Therefore, these 
mode equations asymptotically become. 



dr"^ dr a^(r) 







(8.17) 



The general solutions are given by. 
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/""(r) = exp 



2 



-hr 



-hr 



(8.18) 



where J- 



3/2 1 



and A^' 



3/2 1 



are the Bessel and Neumann functions, respectively. The coeffi- 
cients, c±{q) can be computed for large q by matching fq^^ir) with the WKB approximation 
to the exact mode functions /^(t) that obey the initial conditions (|5.19| ). The WKB ap- 
proximation to fqir) has been computed in ref. [Q, and we find for large q, 



( \ 



l--{h + A) + 0{q-^) 



nh 






/ 8^ 


H 











where 



Jo 



(8.19) 



(8.20) 



In the r ^ cxD limit, we have for fixed q, 




c_{q) 



(8.21) 



which are independent of time asymptotically, and explains why the power spectrum of 
quantum fluctuations freezes at times larger than the spinodal. This behavior is conffimed 
numerically: fig. 14 shows In |/g(T)p vs. r for g = 0,4, 10. Clearly at early times the mode 
functions grow exponentially, and at times of the order of r^, when (7S(r) ~ 1 the mode 
functions freeze-out and become independent of time. Notice that the largest q modes have 
grown the least, explaining why the integral is dominated by g < 10 — 20. 
For asymptotically large times, g^ir) is given by. 



gT.(oo) = g [ — coth 

Jo Q 



2T 



2h 



IT 



(8.22) 



where only one term in the UV subtraction survived in the r = oo limit. The factor coth 
in eg. ( p. 22 ) takes into account the nonzero initial temperature T. 

For consistency, this integral must converge and be equal to 1 as given by the sum rule. 
For this to be the case and to avoid the potential infrared divergence in (|8.22|) , the coefficients 
c_(g) must vanish at g = 0. The mode functions are finite in the g ^ limit provided. 



c„(g) ''=°C g 



3/2 



(8.23) 



where C is a constant. 

The numerical analysis and fig. 12 clearly show that the mode functions remain finite as 
g — > 0, and the coefficient C can be read off from these figures. This is a remarkable result. 
It is well known that for free massless minimally coupled fields in de Sitter space-time with 
Bunch-Davies boundary conditions, the fluctuation contribution {ilj'^{x,t)) grows linearly 
in time as a consequence of the logarithmic divergence in the integrals [^. However, in 
our case, although the asymptotic mode functions are free, the coefficients that multiply 
the Bessel functions of order 3/2 have all the information of the interaction and initial 
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conditions and must lead to the consistency of the sum rule. Clearly the sum rule and the 
initial conditions for the mode functions prevent the coefficients c±{q) from describing the 
Bunch-Davies vacuum. These coefficients are completely determined by the initial conditions 
and the dynamics. This is the reason why the fluctuation freezes at long times unlike in the 
free case in which they grow linearly ||35|| . 

It is easy to see from eqs. (|5.22| )- (|5.23| ) and ( ^.21|) that the energy and pressure vanish for 



r ^ oo. 

Analogously, the two point correlation function can be computed in the late time regime 
using the asymptotic results obtained above. Inserting eq.( 8.1^ ) for the mode functions in 
eq. (|5.26|) yields the asymptotic behavior: 

S{p,T) ^=°° -i- Hqdq sin(gp) coth ( ^) ^ | c„(g) P • (8.24) 



Att'^P Jo \2TJ Ti 

The asymptotic behavior in time of the equal time correlation function is thus solely a 
function of r. The large r behavior of 5'(r, +00) is determined by the singularities of | c_- {k) p 
in the complex k plane. We find an exponential decrease, 

5(p, +00) C '—^ , (8.25) 
P 

where p = i/^ is the pole nearest to the real axis and C is some constant. Thus we see that 
the freeze-out of the mode functions leads to the freeze-out of the correlation length ^. The 
result of the numerical analysis is shown in fig. 16 which confirms this behavior and provides 
the asymptotic value for ^ 1/h. From these figures it is also clear that the freeze-out time 
is given by the expansion time scale, 1/h. More precisely, the numerical values for ^ can be 
accurately reproduced by the following formula obtained by a numerical fit 

~ 1.02 + 0.2 \nh + OMh + ... . 

This situation must be contrasted with that in Minkowski space-time ||18| where the 



correlation length grows as ^{t) ^ ^/T during the stage of spinodal growth. Eventually, this 
correlation length saturates to a fairly large value that is typically several times larger than 
the zero temperature correlation length [IT3. We see that in the de Sitter case the domains 
are always horizon-sized. 



D. Inflaton rolling down (7?(0) 7^ 0): classical or quantum behavior? 

Above we have analyzed the situation when 77(0) = (or in dimensionful variables 
0(0) = 0). The typical analysis of infiaton dynamics in the literature involves the classical 
evolution of 0(t) with an initial condition in which 0(0) is very close to zero (i.e. the top 
of the potential hill) in the 'slow-roll' regime, for which -C 3if0. Thus, it is important 
to quantify the initial conditions on 0(t) for which the dynamics will be determined by 
the classical evolution of 0(t) and those for which the quantum fiuctuations dominate the 
dynamics. We can provide a criterion to separate classical from quantum dynamics by 
analyzing the relevant time scales, estimated by neglecting non-linearities and backreaction 
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effects. We consider the evolution of the zero mode in terms of dimensionless variables, and 
choose 77(0) 7^ and 77(0) = 0. (?)(0) 7^ simply corresponds to a shift in origin of time). 
We assume ?7(0)^ << 1 which is the relevant case where spinodal instabilities are important. 
We find 

ri{r) ^ ri{0) e^'^-l^^^^^ . (8.26) 

The non-linearities will become important and eventually terminate inflation when 
ri{T) ^ 1. This corresponds to a time scale given by 



Comparing this time scale to the spinodal time scale given by (|8.10|) , for which quantum 
fluctuations grow to be of C(l), we see that when, 

r]{0) « g^/\ (8.28) 

the quantum fluctuations will grow to be 0{1) much earlier than the zero mode for Tj > 
(for Tj = the bound becomes 77(0) << g^^"^). In this case the dynamics will be driven 
completely by the quantum fluctuations, as the zero mode will be rolling down the potential 
hill very slowly and will not grow enough to compete with the quantum fluctuations before 
the fluctuations grow to overcome the tree level terms in the equations of motion. In this 
case, as argued previously, the large N and Hartree approximations will be completely 
equivalent during the time scales of interest. 
On the other hand, if 

r]{0) » g'/\ (8.29) 

then the zero mode will roll and become 0{1) before the fluctuations have enough time to 
grow to C(l) (?7(0) >> g^^'^ for Tj = ). In this case, the dynamics will be dominated 
by the rolling of the zero mode and is mostly classical. The quantum fluctuations remain 
perturbatively small throughout the inflationary stage which will end when the velocity of 
the zero mode modifies the equation of state to terminate de Sitter expansion. 

For r]{0) ^ g^^^ (or r]{0) «i g^^'^ for Tj = 0), both the rolling of the zero mode and the 
quantum fluctuations will give contributions of the same order to the dynamics. In this 
case, the quantum fluctuations will be large for the long-wavelength modes and the classical 
approximation to the inflationary dynamics will not be accurate. 

Since the scenario in which t]{0) >> g^^^, in which the dynamics is basically driven by the 
classical evolution of the zero mode has received a great deal of attention in the literature, 
we will not focus on this case, but instead analyze numerically the cases in which //(O) 7^ 
but such that r]{0) < g^^'^. 



E. Numerical Analysis: 



We have evolved the set of equations of motion given by ( ^.18|) numerically with initial 
conditions (|5.19| ) for the large N case, and ( [4. 271) and ([4.28|) , with the corresponding initial 
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conditions ( [4.30| ) on the mode functions for the Hartree case. The numerical code is the 
same as in the previous section with the same relative errors. 
Large N case: 

Figs. 17-18 show ?7(r) and 5'S(r) vs. r for the values g = 10~^^; Tj/Tc = 2; ri{0) = 
10^^; ri{0) = 0. Clearly the dynamics is dominated by the fluctuations; the zero mode 
grows but is always negligible small compared to gJ^i^r). The time scale at which 5'S(r) 
grows to be of order one is about the same as in the case, ri{0) = 0, and all the behavior for 
the mode functions, correlation length, energy density, pressure, etc. is similar to the case 
analyzed in the previous section. 

Asymptotically, we find that the sum rule ( p.20| ) is satisfied to our numerical accuracy. 
This is the same as the situation in Minkowski space-time [|l^,0, and when rj ^ 0, this sum 
rule is nothing but the Ward identity associated with Goldstone's theorem. The fluctuations 
are Goldstone bosons, minimally coupled, and the symmetry is spontaneously broken with 
a very small expectation value for the order parameter as can be read off from fig. 17. For 
T > Tg, the dynamics freezes completely and the zero mode and the fluctuations achieve 
their asymptotic values much in the same way as in the case rj = studied in the previous 
section. Again, the correlation length becomes independent of time with ^{oo) ~ 1//;, in a 
time scale given by 1/h. 

Because there is a damping term in the zero mode equation, it is reasonable to assume 
that asymptotically there will be a solution with a constant value of rj. Then the Ward- 
identity, 77(00) [—1 + 77^(00) + (?S(oo)] = 0, must be fulfilled. In the large case, the only 
stationary solutions are i) r] = 0; gT,{oo) = 1, or, ii) 77(00) 7^ 0; — 1 + 77^(00) + (?S(oo) = 0. 
To have a consistent solution of the mode functions, it must be that the effective mass term 
[— 1 + 77^(00) -|- gT,{oo)] vanishes asymptotically, leading to the mode equations for massless, 
minimally coupled modes which are asymptotically independent of time as shown in the 



previous section (see eg. (|8.21|) ). Furthermore, from fig.18 it is clear that g^ir) remains 



constant at long times, again unlike the case of free massless fields with Bunch-Davies 
boundary conditions in which case the fluctuation grows linearly in time ||35[| . 
Hartree case: 

Figure 17 also shows the evolution of the zero mode in the large A^ and Hartree case. 
Although there is a quantitative difference in the amplitude of the zero mode, in both cases 
it is extremely small and gives a negligible contribution to the dynamics. In the Hartree 
case, however, there is no equivalent of the large A^ sum rule; the only stationary solution 
for ?7 7^ is, 77^(00) = 1; gT,{oo) = 0. Such a solution leads to mode equations with a 
positive mass term and mode functions that vanish exponentially fast for r ^ 00 for all 
momenta. However, whether the asymptotic behavior of the Hartree solution is achieved 
within the interesting time scales is a matter of initial conditions. For example in fig. 17 
the initial condition is such that the time scale for growth of the quantum fluctuations is 
much shorter than the time scale for which the amplitude of the zero mode grows large 
and the non-linearities become important. In the large A^ case the sum rule is satisfied 
with a large value of the quantum fluctuations. In the Hartree case the equivalent sum rule 
— 1 + STjfj + 3gT,H = is satisfied for a very small rin and a gJ^H ~ 1/3. The modes become 
effectively massless and they stop growing. 
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FIG. 17. r/(r) vs. r for A = IQ-^^. r = 2; h = 2; r?(0) 
(solid curve) and Hartree (dashed curve). 



10"^; ?7(0) = for large N 
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FIG. 18. 5S(t) vs. r for the same values of the parameters as in fig. 17, for large (solid 
curve) and Hartree (dashed curve). 



The equation for the zero mode (see eq.( 4.27|) ) still has an uncancelled piece of the non- 
linearity, — 2?7^; however the derivatives and the amplitude of rj are all extremely small and 
though the zero mode still evolves in time, it does so extremely slowly. In fact the Hartree 
curve in fig. 17 has an extremely small positive slope asymptotically, and while tjh grows very 
slowly, q'Eh diminishes at the same rate. In the case shown in fig. 17, we find numerically 
that tih/vh ~ 10^'' at r = 150. Before this time most of the interesting dynamics that can 
be captured with a fixed de Sitter background had already taken place, and the backreaction 
of the fluctuations on the metric becomes substantial requiring an analysis that treats the 
scale factor dynamically. 

The conclusion of our analysis is that in the region of initial conditions for which the 
quantum fluctuations dominate the dynamics, that is for r]{0) << (?^/^, both large and 
Hartree give the same answer on the relevant time scales. The figures for 7i{T)/7i{0) are 
numerically indistinguishable from the case of figs. 1. 

We see that in the large case the zero mode rolls to a final amplitude which is 0{1) 
and of the same order as gT,{oo) and the sum rule is satisfied. However, the Hartree case 
clearly shows the asymptotics analyzed above with rjn (oo) = l;gT,H{oo) = 0. 

This particular borderline case is certainly not generic and would imply some fine tuning 
of initial conditions. Finally the case in which 7](0) >> g^^^ (or g^^"^ for Tj = 0) is basically 
classical in that the dynamics is completely given by the classical rolling of the zero mode 
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and the fluctuations are always perturbatively smalL 



F. Discussion and Conclusions for the de Sitter background 

We have identifled analytically and numerically two distinct regimes for the dynamics 
determined by the initial condition on the expectation value of the zero mode of the inflaton 



m 



1. When 77(0) << g^^^ (or g^^"^ for Tj = 0), the dynamics is driven by quantum (and 
thermal) fluctuations. Spinodal instabilities grow and eventually compete with tree 
level terms at a time scale, > —3hln[g]/2. The growth of spinodal fluctuations 
translates into the growth of spatially correlated domains which attain a maximum 
correlation length (domain size) of the order of the horizon. For very weak coupling 
and h > 1 this time scale can easily accommodate enough e-folds for inflation to solve 
the flatness and horizon problems. The quantum fluctuations modify the equation 
of state dramatically providing a means for a graceful exit to the inflationary stage 
without slow-roll. 

This non-perturbative description of the non-equilibrium effects in this regime in which 
quantum (and thermal) fluctuations are most important is borne out by both the large 
and Hartree approximations. Thus our analysis provides a reliable understanding 
of the relevant non-perturbative, non-equilibrium effects of the fluctuations that have 
not been revealed before in this setting - |T6[ . 



These initial conditions are rather natural if the de Sitter era arises during a phase 
transition from a radiation dominated high temperature phase in local thermodynamic 
equilibrium, in which the order parameter and its time derivative vanish. 

When 77(0) >> g^^^ (or g^^"^ for Tj = 0), the dynamics is driven solely by the classical 
evolution of the inflaton zero mode. The quantum and thermal fluctuations are always 
perturbatively small (after renormalization), and their contribution to the dynamics is 
negligible for weak couplings. The de Sitter era will end when the kinetic contribution 
to the energy becomes of the same order as the 'vacuum' term. This is the realm 
of the slow-roll analysis whose characteristics and consequences have been analyzed 
in the literature at length. These initial conditions, however, necessarily imply some 
initial state either with a biasing fleld that favors a non-zero initial expectation value, 
or that in the radiation dominated stage, prior to the phase transition, the state was 
strongly out of equilibrium with an expectation value of the zero mode different from 
zero. Although such a state cannot be ruled out and would naturally arise in chaotic 
scenarios, the description of the phase transition in this case requires further input on 
the nature of the state prior to the phase transition. 

IX. SELF-CONSISTENT EVOLUTION OF MATTER FIELDS WITH A 
DYNAMICAL COSMOLOGICAL BACKGROUND 



We present in this section the full self-consistent matter-geometry dynamics [0]. That 



is, the scale factor a(t) is here a dynamical variable determined by the Einstein-Friedman 
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eq.(|]2T])-(|5]2|) coupled with the matter evolution eqs.(|]T8D-(|]T9]). 



In order to provide the full solution we now must provide the values of ^^(O), ?7(0), and 
hf). Assuming that the inflationary epoch is associated with a phase transition at the GUT 
scale, this requires that Nm'^/g ^ (10^^ Gev )^ and assuming the bound on the scalar self- 
coupling g ^ 10~^^ — 10~^^ (this will be seen later to be a compatible requirement), we find 
that ho ^ N'^^'^ which we will take to be reasonably given by /iq ~ 1 — 10 (for example in 
popular cut's ^ 20 depending on particular representations). 

We will begin by studying the case of most interest from the point of view of describing 
the phase transition: //(O) = and 77(0) = 0, which are the initial conditions that led to 
puzzling questions. With these initial conditions, the evolution equation for the zero mode 
eq. (|5.18|) determines that //(r) = by symmetry. 



A. Early time dynamics: 

Before engaging in the numerical study, it proves illuminating to obtain an estimate of 
the relevant time scales and an intuitive idea of the main features of the dynamics. Because 
the coupling is so weak {g ~ 10^^^ ^ 1) and after renormalization the contribution from 
the quantum fluctuations to the equations of motion is finite, we can neglect all the terms 
proportional to g in eqs.( ^.22D and ( p. 181 ). 

For the case where we choose ri{T) = and the evolution equations for the mode functions 
are those for an inverted oscillator in De Sitter space-time, which have been studied in sec. 



VIII One obtains the approximate solutions ( ^.3|) -( 

After the physical wavevectors cross the horizon, i.e. when qe~^"'^/ho ^ 1 we find that 
the mode functions factorize: 

This result reveals a very important feature: because of the negative mass squared term 
in the matter Lagrangian leading to symmetry breaking (and u > 3/2), we see that all 
of the mode functions grow exponentially after horizon crossing (for positive mass squared 
u < 3/2, and they would decrease exponentially after horizon crossing). This exponential 
growth is a consequence of the spinodal instabilities which is a hallmark of the process of 
phase separation that occurs to complete the phase transition. We note, in addition that 
the time dependence is exactly given by that of the q = mode, i.e. the zero mode, which is 
a consequence of the redshifting of the wavevectors and the fact that after horizon crossing 
the contribution of the term q'^/a^{T) in the equations of motion become negligible. We 
clearly see that the quantum fluctuations grow exponentially and they will begin to be of 
the order of the tree level terms in the equations of motion when 5'S(r) ^ 1. At large times 

S(r) ^ ^2(^^)^2g(2.-3)ho. ^ 

with J-'{ho) a finite constant that depends on the initial conditions and is found numerically 



to be of C(l) [see fig.|4|. 
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In terms of the initial dimensionful variables, the condition 5'S(r) ^ 1 translates to 
< ip'^{x,t) >H~ 2m^/5f, i.e. the quantum fluctuations sample the minima of the (renormal- 
ized) tree level potential. We find that the time at which the contribution of the quantum 
fluctuations becomes of the same order as the tree level terms is estimated to be 



1 



(2i/ - 3)ho 



In 



1 



ghlJ^^{ho 



-hn In 



ghlJ^^ho) 



(9.2) 



At this time, the contribution of the quantum fluctuations makes the back reaction very im- 
portant and, as will be seen numerically, this translates into the fact that also determines 
the end of the De Sitter era and the end of inflation. The total number of e-folds during the 
stage of exponential expansion of the scale factor (constant ho) is given by 



2z/-3 



In 



g hi J^^ho) 



hi In 



q hi J^^iho) 



+ 0{1) 



(9.3) 



For large /iq we see that the number of e-folds scales 8bS clS well as with the logarithm of 



the inverse coupling. These results (|9.1| - |9.3| ) will be confirmed numerically below and will be 
of paramount importance for the interpretation of the main consequences of the dynamical 
evolution. 

As discussed in sec. VIII. D, the early time dynamics is dominated by classical or quantum 
effects depending on the ratio between the time scales Tc and r^. 

If Tc is much smaller than the spinodal time Tg given by eq.( |9.2|) then the classical 
evolution of the zero mode will dominate the dynamics and the quantum fluctuations will 
not become very large, although they will still undergo spinodal growth. On the other hand, 
if Tc 2> Tg the quantum fluctuations will grow to be very large well before the zero mode 
reaches the non-linear regime. In this case the dynamics will be determined completely by 
the quantum fluctuations. Then the criterion for the classical or quantum dynamics is given 
by 

77(0) ^ y/g ho =^ classical dynamics 

77(0) ^ ho =^ quantum dynamics (9.4) 

or in terms of dimensionful variables 0(0) ^ Hq leads to classical dynamics and 0(0) <^ Hq 
leads to quantum dynamics. 

However, even when the classical evolution of the zero mode dominates the dynamics, 
the quantum fluctuations grow exponentially after horizon crossing unless the value of 0(t) 
is very close to the minimum of the tree level potential. In the large N approximation the 
spinodal line, that is the values of 0(t) for which there are spinodal instabilities, reaches all 
the way to the minimum of the tree level potential as can be seen from the equations of 
motion for the mode functions. Therefore even in the classical case one must understand 
how to deal with quantum fluctuations that grow after horizon crossing. 



B. Numerics 

The time evolution is carried out by means of a fourth order Runge-Kutta routine with 
adaptive step-sizing while the momentum integrals are carried out using an 11-point Newton- 
Cotes integrator. The relative errors in both the differential equation and the integration 
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are of order 10~^. We find that the energy is covariantly conserved throughout the evolution 
to better than a part in a thousand. Figs. |l9|-0 show (y'S(r) vs. r, /i(r) vs. r and In |/g(r)p 
vs. r for several values of q with larger q' s corresponding to successively lower curves. Figs. 
22| , p3| show p{T)/e{T) and the horizon size h~^{T) for g = 10~^^ ; ri{0) = ; 77(0) = and 
we have chosen the representative value ho = 2.0. 



1 .0 - 

^ 0.5 - 



FIG. 19. c/S(t) vs. r, for 7/(0) = 0, 77(0) = 0,g = IQ-^'^, ho = 2.0 



FIG. 20. H{t) vs. t, for 77(0) = 0,7)(0) =0,g = lO-^^^o = 2.0 




FIG. 21. ln|/g(T)r vs. r, for 7/(0) = 0,7/(0) = 0,g = 10-l^/^o = 2.0 for q = 0.0,5,10,15,20 
with smaller q corresponding to larger values of In |/g(T)p. 



Figs. |T9| and ^ show clearly that when the contribution of the quantum fluctuations 
(/S(r) becomes of order 1 inflation ends, and the time scale for (/S(r) to reach 0{1) is very 
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well described by the estimate ( |9.2|) . From fig. 19 we see that this happens for t = Ts 90, 
leading to a number of e-folds ~ 180 which is correctly estimated by (|9^ , [97^ ). 

Fig. |21] shows clearly the factorization of the modes after they cross the horizon as 
described by eg. ( |9.1|) . The slopes of all the curves after they become straight lines in fig.0 is 



given exactly by (2z/ — 3), whereas the intercept depends on the initial condition on the mode 
function and the larger the value of q the smaller the intercept because the amplitude of the 
mode function is smaller initially. Although the intercept depends on the initial conditions 
on the long- wavelength modes, the slope is independent of the value of q and is the same as 
what would be obtained in the linear approximation for the square of the zero mode at times 
long enough that the decaying solution can be neglected but short enough that the effect 
of the non-linearities is very small. Notice from the figure that when inflation ends and the 
non-linearities become important all of the modes effectively saturate. This is also what one 
would expect from the solution of the zero mode: exponential growth in early-intermediate 
times (neglecting the decaying solution), with a growth exponent given by (z/ — 3/2) and an 
asymptotic behavior of small oscillations around the equilibrium position, which for the zero 
mode is ?7 = 1, but for the q ^ ^ modes depends on the initial conditions. All of the mode 
functions have this behavior once they cross the horizon. We have also studied the phases 
of the mode functions and we found that they freeze after horizon crossing in the sense that 
they become independent of time. This is natural since both the real and imaginary parts 
of /^(t) obey the same equation but with different boundary conditions. After the physical 
wavelength crosses the horizon, the dynamics is insensitive to the value of q for real and 
imaginary parts and the phases become independent of time. Again, this is a consequence 
of the factorization of the modes. 

The growth of the quantum fluctuations is sufficient to end inflation at a time given by 
Ts in eq.( |9.2| ). Furthermore fig. shows that during the inflationary epoch p{t)/6{t) ^ —1 
and the end of inflation is rather sharp at with p(r)/£:(r) oscillating between ±1 with 
zero average over the cycles, resulting in matter domination. Fig. 12^ shows this feature very 
clearly; h{T) is constant during the de Sitter epoch and becomes matter dominated after 
the end of inflation with h~^{T) ^ |(r — r^). There are small oscillations around this value 
because both p(r) and e(r) oscillate. These oscillations are a result of small oscillations 
of the mode functions after they saturate, and are also a feature of the solution for a zero 
mode. 

All of these features hold for a variety of initial conditions. As an example, we show in 
ref. W^ the case of an initial Hubble parameter of ho = 10. 
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C. Zero Mode Assembly: 

This remarkable feature of factorization of the mode functions after horizon crossing can 
be elegantly summarized as 

/fc(^)lfcp^t)«// = aiq, h)foiT), (9.5) 

with kphit) = ke^^^ being the physical momentum, g[q, h) a complex constant, and /o(t) a 
real function of time that satisfies the mode equation with g = and real initial conditions 
which will be inferred later. Since the factor g{q, h) depends solely on the initial conditions 
on the mode functions, it turns out that for two mode functions corresponding to momenta 
ki, k2 that have crossed the horizon at times ti > t2, the ratio of the two mode functions at 
time t, {ts > t > ti > ^2) is 

Then if we consider the contribution of these modes to the renormalized quantum fiuctuations 
a long time after the beginning of inflation (so as to neglect the decaying solutions), we flnd 
that 

c/E(t) ^ Ce^'^^-^'^h^ + small , 
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where 'small' stands for the contribution of mode functions associated with momenta that 
have not yet crossed the horizon at time r, which give a perturbatively small (of order g) 
contribution. We find that several e-folds after the beginning of inflation but well before 
inflation ends, this factorization of superhorizon modes implies the following: 



gJq'dq\f^{T)\^\Co\'f^iT), 



9 



9 I q'dq\f'^ir)\^\Co\'f'oir), 
2/ ^dq |/2(r) 



-foir), 



(9.6) 
(9.7) 

(9.8) 



where we have neglected the weak time dependence arising from the perturbatively small 
contributions of the short- wavelength modes that have not yet crossed the horizon, and the 
integrals above are to be understood as the fully renormalized (subtracted), finite integrals. 
For ?7 = 0, we note that (|9l6|) and the fact that /o(t) obeys the equation of motion for the 
mode with g = leads at once to the conclusion that in this regime [^'^(r)]^ = |Co|/o(t) 
obeys the zero mode equation of motion 



\Co\foir) = 



(9.9) 



It is clear that several e-folds after the beginning of inflation, we can define an effective 
zero mode as 



(7S(r), or in dimensionful variables, 0e//(^) = {'4''^{x,t)) r 



(9.10) 



Although this identification seems natural, we emphasize that it is by no means a trivial or 
ad- hoc statement. There are several important features that allow an unambiguous identifi- 
cation: i) [{tp'^ {x , t)) n] is a fully renormalized operator product and hence finite, ii) because 
of the factorization of the superhorizon modes that enter in the evaluation of [{tp'^{x,t)) r], 
^effif) (|9.10|) obeys the equation of motion for the zero mode, iii) this identification is valid 



several e-folds after the beginning of inflation, after the transient decaying solutions have 
died away and the integral in {ip'^{x,t)) is dominated by the modes with wavevector k that 
have crossed the horizon at t{k) <^ t. Numerically we see that this identification holds 
throughout the dynamics except for a very few e-folds at the beginning of infiation. This 
factorization determines at once the initial conditions of the effective zero mode that can be 
extracted numerically: after the first few e-folds and long before the end of infiation we find 



(t)eff{t) = (peffiO) e 



(9.11) 



where we parameterized 



to make contact with the literature. As is shown in fig. (|24D , we find numerically that 
J-'{H/m) ~ for a large range of 0.1 < H/m < 50 and that this quantity depends on 
the initial conditions of the long wavelength modes. 
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h 



FIG. 24. J^{H/m) vs. H, where J^{H/m) is defined by the relation 0e//(O) = {H /2tt)F{H /m) 
(see eqs. ( plO| ) and (CT))- 



h.,(t) 




FIG. 25. r?2^/r) vs. r, for r?e//(0) = 3.94 x 10-^7)e//(0) = 0.317r/e//(0), 5 = 10-^^ /iq = 2.0. 
The initial conditions were obtained by fitting the intermediate time regime of (7S(r) in fig. 19 
T]eff{T) is the solution of eq.( 9.12| ) with these initial conditions. 



FIG. 26. /i(r) vs. r, obtained from the solution of eqs. ( 9.12| ) and ( |9.15| ) with the conditions 



of fig. 25 



Therefore, in summary, the effective composite zero mode obeys 



d' 



d 



rieff{r) = ; r]eff{r = 0) = {u - -) r/e//(0) , 



(9.12) 



where i]eff{0) = ^^^^ 0e//(O) is obtained numerically for a given by fitting the interme- 



60 



diate time behavior of 5'S(r) with the growing zero mode solution. Recall that Xr = Sir'^ g. 

Furthermore, this analysis shows that in the case rj = 0, the renormalized energy and 
pressure in this regime in which the renormalized integrals are dominated by the superhori- 
zon modes are given by 

^«(^) - ^ ^eff + i + Vhf) I ' (9-13) 

{p + e)n^'-^{v!ff}, (9.14) 



R 



where we have neglected the contribution proportional to l/a^(r) because it is effectively 
red-shifted away after just a few e-folds. We found numerically that this term is negligible 
after the interval of time necessary for the superhorizon modes to dominate the contribution 
to the integrals. Then the dynamics of the scale factor is given by 
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We have numerically evolved the set of effective equations ( |9.12| , |9.15| ) by extracting the 



initial condition for the effective zero mode from the intermediate time behavior of gT,(T). 
We found a remarkable agreement between the evolution of r^g^j and 5'S(r) and between the 
dynamics of the scale factor in terms of the evolution of rjeffi^r), and the full dynamics of 
the scale factor and quantum fluctuations within our numerical accuracy. Figs. and p6| 
show the evolution of rjlffir) and hi^r) respectively from the classical evolution eqs. (|9.12|) 
and (|9.15| ) using the initial condition i]eff{0) extracted from the exponential fit of 5'S(r) in 



the intermediate regime. These figures should be compared to figs. |T9| and We have also 
numerically compared p/e given solely by the dynamics of the effective zero mode and it is 
again numerically indistinguishable from that obtained with the full evolution of the mode 
functions. 

This is one of the main results of our work. In summary: the modes that become 
superhorizon sized and grow through the spinodal instabilities assemble themselves into an 
effective composite zero mode a few e-folds after the beginning of inflation. This effective 
zero mode drives the dynamics of the FRW scale factor, terminating inflation when the non- 
linearities become important. In terms of the underlying fluctuations, the spinodal growth of 
superhorizon modes gives a non-perturbatively large contribution to the energy momentum 
tensor that drives the dynamics of the scale factor. Inflation terminates when the mean 
square root fluctuation probes the equilibrium minima of the tree level potential. 

This phenomenon of zero mode assembly, i.e. the 'classicalization' of quantum mechan- 
ical fluctuations that grow after horizon crossing is very similar to the interpretation of 
'decoherence without decoherence' of Starobinsky and Polarski p3|| . 



The extension of this analysis to the case for which 77(0) 7^ is straightforward. Since 



both ■r]{T) and ^ g'Z{T) = |Co|/o(t) obey the equation for the zero mode, eq.( ^.18|) , it is clear 
that we can generalize our definition of the effective zero mode to be 

Veffir)^^r^^iT)+gJ:iT) . (9.16) 
which obeys the equation of motion of a classical zero mode: 
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Veffir) = . 



(9.17) 



If this effective zero mode is to drive the FRW expansion, then the additional condition 



(9.18) 



must also be satisfied. One can easily show that this relation is indeed satisfied if the 
mode functions factorize as in (|9.5|) and if the integrals ( |9.6D - ( |9.8|) are dominated by the 
contributions of the superhorizon mode functions. This leads to the conclusion that the 
gravitational dynamics is given by eqs. (|9.13| ) - ( |9.15|) with rjeffij) defined by (|9.16| ). 

We see that in all cases, the full large N quantum dynamics in these models of inflation- 
ary phase transitions is well approximated by the equivalent dynamics of a homogeneous, 
classical scalar field with initial conditions on the effective field ?7e//(0) > ^hoJ^{ho). We 
have verified these results numerically for the field and scale factor dynamics, finding that 
the effective classical dynamics reproduces the results of the full dynamics to within our 
numerical accuracy. We have also checked numerically that the estimate for the classical to 
quantum crossover given by eg. (|9.4| ) is quantitatively correct. Thus in the classical case in 
which ri{0) ^ ^ ho we find that Tjeff^r) = rj^r), whereas in the opposite, quantum case 



Veffir) = ^^S(r). 

This remarkable feature of zero mode assembly of long-wavelength, spinodally unstable 
modes is a consequence of the presence of the horizon. It also explains why, despite the fact 
that asymptotically the fluctuations sample the broken symmetry state, the equation of state 
is that of matter. Since the excitations in the broken symmetry state are massless Goldstone 
bosons one would expect radiation domination. However, the assembly phenomenon, i.e. 
the redshifting of the wave vectors, makes these modes behave exactly like zero momentum 
modes that give an equation of state of matter (upon averaging over the small oscillations 
around the minimum). 

Subhorizon modes at the end of inflation with q > e'^"^' do not participate in the 
zero mode assembly. The behavior of such modes do depend on q after the end of inflation. 
Notice that these modes have extremely large comoving q since ho e'*"'^" > 10^^. As discussed 



in sec. VII such modes decrease with time after inflation as ~ l/a{T) |14|. 



D. Making sense of small fluctuations: 

Having recognized the effective classical variable that can be interpreted as the compo- 
nent of the field that drives the FRW background and rolls down the classical potential hill, 
we want to recognize unambiguously the small fluctuations. We have argued above that 
after horizon crossing, all of the mode functions evolve proportionally to the zero mode, and 
the question arises: which modes are assembled into the effective zero mode whose dynamics 
drives the evolution of the FRW scale factor and which modes are treated as perturbations? 
In principle every k ^ mode provides some spatial inhomogeneity, and assembling these 
into an effective homogeneous zero mode seems in principle to do away with the very in- 
homogeneities that one wants to study. However, scales of cosmological importance today 
first crossed the horizon during the last 60 or so e-folds of inflation. Recently Grishchuk 
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4^ has argued that the sensitivity of the measurements of AT/T probe inhomogeneities 



on scales ~ 500 times the size of the present horizon. Therefore scales that are larger than 
these and that have first crossed the horizon much earlier than the last 60 e-folds of inflation 
are unobservable today and can be treated as an effective homogeneous component, whereas 
the scales that can be probed experimentally via the CMB inhomogeneities today must be 
treated separately as part of the inhomogeneous perturbations of the CMB. 

Thus a consistent description of the dynamics in terms of an effective zero mode plus 
'small' quantum fluctuations can be given provided the following requirements are met: 
a) the total number of e-folds ^ 60, b) all the modes that have crossed the horizon 
before the last 60-65 e-folds are assembled into an effective classical zero mode via 4>effit) = 

[00 (^) + r]^ y c) the modes that cross the horizon during the last 60-65 e-folds 

are accounted as 'small' perturbations. The reason for the requirement a) is that in the 
separation (p{x,t) = 0e//(^) + S(j){x,t) one requires that 5(t){x,t) / (peffif) -C 1. As argued 
above, after the modes cross the horizon, the ratio of amplitudes of the mode functions 
remains constant and given by e^''~2)^^ with AA^ being the number of e-folds between the 
crossing of the smaller k and the crossing of the larger k. Then for t) to be much smaller 
than the effective zero mode, it must be that the Fourier components of correspond to 
very large /c's at the beginning of inflation, so that the effective zero mode can grow for a 
long time before the components of 50 begin to grow under the spinodal instabilities. In 
fact requirement a) is not very severe; in the figs. (19-23) we have taken = 2.0 which 
is a very moderate value and yet for g = 10~^^ the inflationary stage lasts for well over 
100 e-folds, and as argued above, the larger ho for fixed g, the longer is the inflationary 
stage. Therefore under this set of conditions, the classical dynamics of the effective zero 
mode (peffit) drives the FRW background, whereas the inhomogeneous fluctuations S(j){x, t), 
which are made up of Fourier components with wavelengths that are much smaller than the 
horizon at the beginning of inflation and that cross the horizon during the last 60 e-folds, 
provide the inhomogeneities that seed density perturbations. 



E. Scalar Metric Perturbations: 

Having identified the effective zero mode and the 'small perturbations', we are now in 
position to provide an estimate for the amplitude and spectrum of scalar metric perturba- 
tions. We use the clear formulation in ref. |^ in terms of gauge invariant variables. In 
particular we focus on the dynamics of the Bardeen potential |^^, which in longitudinal 
gauge is identified with the Newtonian potential. The equation of motion for the Fourier 
components (in terms of comoving wavevectors) for this variable in terms of the effective 
zero mode is lisll 



(l>eff{t)_ 



+ 2 H{t) 



$fc = 0. 



(9.19) 



We are interested in determining the dynamics of for those wavevectors that cross 
the horizon during the last 60 e-folds before the end of inflation. During the inflationary 
stage the numerical analysis yields to a very good approximation 



H{t) ^ Ho ; 0e//(t) = 0e//(O) e 



(9.20) 
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where Hq is the value of the Hubble constant during inflation, leading to 



(2) 



V 



(9.21) 



The coeflicients 0^,6^ are determined by the initial conditions. 

Since we are interested in the wavevectors that cross the horizon during the last 60 
e-folds, the consistency for the zero mode assembly and the interpretation of 'small pertur- 
bations' requires that there must be many e-folds before the last 60. We are then considering 
wavevectors that were deep inside the horizon at the onset of inflation. is related to 

the canonical 'velocity fleld' that determines scalar perturbations of the metric and which is 
quantized with Bunch-Davies initial conditions for the large /c-mode functions. The relation 
between and v and the initial conditions on v lead at once to a determination of the 



coefficients and 6^ ioi k >> Hq 



in 



^eff 



(0)i 



TC 



2Hok 







(9.22) 



Thus we find that the amplitude of scalar metric perturbations after horizon crossing is 
given by 



\6k{t)\ = k-2\Mt)\ 



^eff 



(0) 



2Hc 

k 



,i2u-3)Hot 



(9.23) 



The power spectrum per logarithmic k interval is given by |5fc(t)p. The time dependence 
of |5fc(t)| displays the unstable growth associated with the spinodal instabilities of super- 
horizon modes and is a hallmark of the phase transition. This time dependence can be also 
understood from the constraint equation that relates the Bardeen potential to the gauge 
invariant fleld fluctuations which in longitudinal gauge are identifled with 6(j){x,t). 
The constraint equation and the evolution equations for the gauge invariant scalar fleld 
fluctuations are [EH 
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+ 3H- + — +M^ 



dt 



4 0e// <i>fe + 2V'{(j)eff) = 



(9.25) 



Since the right hand side of (|9.24|) is proportional to (f)eff/Mpi <^ 1 during the inflationary 
epoch in this model, we can neglect the terms proportional to and $fc on the left hand 
side of (|9.25|) , in which case the equation for the gauge invariant scalar fleld fluctuation is 
the same as for the mode functions. In fact, since 50^* is gauge invariant we can evaluate 
it in the longitudinal gauge wherein it is identifled with the mode functions fk{t). Then 
absorbing a constant of integration in the initial conditions for the Bardeen variable, we flnd 



to 



ait') 0e//(O/fc(O dt' + O 



Ml 



(9.26) 
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and using that oc e^*^"^/^)^"* and that after horizon crossing fk(t) oc e^'^~^/'^^^°^ , one 
obtains at once the time dependence of the Bardeen variable after horizon crossing. In 
particular the time dependence is found to be oc ei'^>^-^)^ot _ j^. ^]^g^ clear that the time 
dependence is a reflection of the spinodal (unstable) growth of the superhorizon field fluc- 
tuations. 

To obtain the amplitude and spectrum of density perturbations at seconc? horizon crossing 



we use the conservation law associated with the gauge invariant variable |^| 



k — TT-T- ; r "Pa 



3 1+p/e 



6 







(9.27) 



which is valid after horizon crossing of the mode with wavevector k. Although this conserva- 
tion law is an exact statement of superhorizon mode solutions of eq.( |9.19D , we have obtained 
solutions assuming that during the inflationary stage H is constant and have neglected the 
H term in Eq. (9.iy). Since during the inflationary stage, 



m = 



4:71 

'Mil 



(9.28) 



and ~ Hq, the above approximation is justified. We then see that (pljrj^t) oc e^'^'^~^^^°^ 
which is the same time dependence as that of ^k{t)- Thus the term proportional to 1/(1 + 
p/e) in Eq. (|9.27|) is indeed constant in time after horizon crossing. On the other hand, 
the term that does not have this denominator evolves in time but is of order (1 + p/e) = 
—2H/3H^ <^ 1 with respect to the constant term and therefore can be neglected. Thus, we 
confirm that the variable ^ is conserved up to the small term proportional to (1 +p/e)^k 
which is negligible during the inflationary stage. This small time dependence is consistent 
with the fact that we neglected the H term in the equation of motion for $fc(t). The validity 
of the conservation law has been recently studied and conflrmed in different contexts |^ . 
Notice that we do not have to assume that $fc vanishes, which in fact does not occur. 

However, upon second horizon crossing it is straightforward to see that $fc(t/) ~ 0. The 
reason for this assertion can be seen as follows: eq.( p.25| ) shows that at long times, when 
the effective zero mode is oscillating around the minimum of the potential with a very small 
amplitude and when the time dependence of the fluctuations has saturated (see flg.3), $fc 
will redshift as ~ l/a{t) and its derivative becomes extremely small. 

Using this conservation law, assuming matter domination at second horizon crossing, 
and $fc(i/) ~ we flnd 



Mtf)\ 



12r(z/)v^ f2H^ 



5 (i^ - f ) J^iHo/m) V k 



(9.29) 



where J-'{Ho/m) determines the initial amplitude of the effective zero mode ( |9.11|) . We can 
now read the power spectrum per logarithmic k interval 



Vsik) = \Sk\' oc A;-2('^-t). 



(9.30) 



leading to the index for scalar density perturbations 
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n. = l-2(.-|) 



(9M) 



For Hq/itl ^ 1, we can expand z/ — 3/2 as a series in m? /Hq in eg .( p .291) . Given that 
the comoving wavenumber of the mode which crosses the horizon n e-folds before the end of 
inflation is = i^oe'-^''""'' where A'^e is given by (|9.3| ), we arrive at the following expression 
for the amplitude of fluctuations on the scale corresponding to n in terms of the De Sitter 
Hubble constant and the coupling A = Svr^^f: 



\Utf)\ 



5V2m^ 



(9.32) 



Here, 7 is Euler's constant. Note the explicit dependence of the amplitude of density per- 
turbations on y/g. For n ^ 60, the factor exp(nm^/3ifo) is (9(100) for Hq/ui = 2, while it 
is 0{1) for Hq/tti > 4. Notice that for Hq/tti large, the amplitude increases approximately 
as {Ho/m)^, which will place strong restrictions on g in such models. 

We remark that we have not included the small corrections to the dynamics of the 
effective zero mode and the scale factor arising from the non-linearities. We have found 
numerically that these nonlinearities are only significant for the modes that cross about 60 
e-folds before the end of inflation for values of the Hubble parameter HQ/mji > 5. The effect 
of these non-linearities in the large N limit is to slow somewhat the exponential growth of 
these modes, with the result of shifting the power spectrum closer to an exact Harrison- 
Zeldovich spectrum with Ug = 1. Since for Ho/rrifi > 5 the power spectrum given by (p.31|) 
differs from one by at most a few percent, the effects of the non-linearities are expected to 
be observationally unimportant. The spectrum given by ( |9.29| ) is similar to that obtained 
in references although the amplitude differs from that obtained there. In addition, we 

do not assume slow roll for which (z/ — |) ^ 1, although this would be the case if ^ 60. 

We emphasize an important feature of the spectrum: it has more power at long wave- 
lengths because u — 3/2 > 0. This is recognized to be a consequence of the spinodal insta- 
bilities that result in the growth of long wavelength modes and therefore in more power for 
these modes. This seems to be a robust prediction of new inflationary scenarios in which the 
potential has negative second derivative in the region of field space that produces inflation. 

It is at this stage that we recognize the consistency of our approach for separating the 
composite effective zero mode from the small fluctuations. We have argued above that many 
more than 60 e-folds are required for consistency, and that the small fluctuations correspond 
to those modes that cross the horizon during the last 60 e-folds of the inflationary stage. 
For these modes Ho/k = e~^°^*^^^ where t*{k) is the time since the beginning of inflation of 
horizon crossing of the mode with wavevector k. The scale that corresponds to the Hubble 
radius today Aq = 27r/fco is the first to cross during the last 60 or so e-folds before the 
end of inflation. Smaller scales today will correspond to k > ko at the onset of inflation 
since they will cross the flrst horizon later and therefore will reenter earlier. The bound 
on \Sko\ oc AT/T < 10^^ on these scales provides a lower bound on the number of e-folds 
required for these type of models to be consistent: 



AT, > 60 + 



12 



ln{u — 



u — 



V — 



(9.33) 
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where we have written the total number of e-folds as = Hq t*{ko) + 60. This in turn can 
be translated into a bound on the coupling constant using the estimate given by eq. (|9.3|) . 
The four year COBE DMR Sky Map |^ gives n ~ 1.2 ± 0.3 thus providing an upper 



bound on u 

0<z/-^< 0.05 (9.34) 

corresponding to ho > 2.6. We then find that these values of ho and A ~ 10^^^ — 10~^^ 
provide sufficient e-folds to satisfy the constraint for scalar density perturbations. 



F. Tensor Metric Perturbations: 



The scalar field does not couple to the tensor (gravitational wave) modes directly, and 
the tensor perturbations are gauge invariant from the beginning. Their dynamical evolution 
is completely determined by the dynamics of the scale factor [EsLEQI. Having established 



numerically that the infiationary epoch is characterized by H / Hq -C 1 and that scales 
of cosmological interest cross the horizon during the stage in which this approximation is 
excellent, we can just borrow the known result for the power spectrum of gravitational waves 



produced during infiation extrapolated to the matter era p^ , |49 



^Tik) ^ . (9.35) 



Thus the spectrum to this order is scale invariant (Harrison-Zeldovich) with an amplitude 
of the order m7AM|,. Then, for values of m f« lO^^ _ iqu q^^ ^nd A ^ lO'^^ _ iq-u 
one finds that the amplitude is < 10^^'^ which is much smaller than the amplitude of scalar 
density perturbations. As usual the amplification of scalar perturbations is a consequence 
of the equation of state during the infiationary epoch. 



G. Contact with the Reconstruction Program: 

The program of reconstruction of the infiationary potential seeks to establish a rela- 
tionship between features of the infiationary scalar potential and the spectrum of scalar and 
tensor perturbations. This program, in combination with measurements of scalar and tensor 
components either from refined measurements of temperature inhomogeneities of the CMB 
or through galaxy correlation functions will then offer a glimpse of the possible realization of 
the inflation [Q^. Such a reconstruction program is based on the slow roll approximation 
and the spectral index of scalar and tensor perturbations are obtained in a perturbative 



expansion in the slow roll parameters |5C.51| 



4,) ^ ^ . (0.36) 

^(0) = • (9.37) 

Hep 
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We can make contact with the reconstruction program by identifying above with our 0e// 
after the first few e-folds of inflation needed to assemble the effective zero mode from the 
quantum fluctuations. We have numerically established that for the weak scalar coupling 
required for the consistency of these models, the cosmologically interesting scales cross the 
horizon during the epoch in which H ^ Hq ; (peff ~ (z^ — 3/2) Hq (peff V ~ 171%/ X ^ 0e//- 
In this case we find 

v{<Peff) = -{v - \) ; e(0e//) ^ 0{\) « r7(0e//). (9.38) 

With these identifications, and in the notation of |]50| , |5T| the reconstruction program 
predicts the index for scalar density perturbations Us given by 

n,-l = -2(^z/-^)+0(A), (9.39) 

which coincides with the index for the power spectrum per logarithmic interval \5k\^ with 
l^fcl given by eq.( |9.29| ). We must note however that our treatment did not assume slow roll 
for which (z/ — |) ^ 1. Our self-consistent, non-perturbative study of the dynamics plus the 
underlying requirements for the identification of a composite operator acting as an effective 
zero mode, validates the reconstruction program in weakly coupled new inflationary models. 



X. DECOHERENCE: QUANTUM TO CLASSICAL TRANSITION DURING 

INFLATION 

An important aspect of cosmological perturbations is that they are of quantum origin 
but eventually they become classical as they are responsible for the small classical metric 
perturbations. This quantum to classical crossover is associated with a decoherence process 
and has received much attention |^J5^ . 



In this section we study the quantum to classical transition of superhorizon modes for 
the Bardeen variable by relating these to the field mode functions through eq. (|9.26|) and 



analyzing the full time evolution of the density matrix of the matter field. Eq.( |9.26D es- 
tablishes that in the models under consideration the classicality of the Bardeen variable is 
determined by the classicality of the scalar field modes. 

In the situation under consideration, long-wavelength field modes become spinodally 
unstable and grow exponentially after horizon crossing. The factorization ( |9.1| ) results in 
the phases of these modes 'freezing out'. This feature and the growth in amplitude entail 
that these modes become classical. Eq.( |9.26| ) in turn implies that these features also apply 
to the superhorizon modes of the Bardeen potential. 

Therefore we can address the quantum to classical transition of the Bardeen variable 
(gravitational potential) by analyzing the evolution of the density matrix for the matter 
field. 



To make contact with previous work |43,52| we choose to study the evolution of the field 



density matrix in conformal time, although the same features will be displayed in comoving 
time. 

In the large N or Hartree (also in the self-consistent one-loop) approximation, the den- 
sity matrix is Gaussian, and defined by a normalization factor, a complex covariance that 
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determines the diagonal matrix elements, and a real covariance that determines the mixing 
in the Schrodinger representation as discussed in ref. |jl3| (and references therein). 
That is, the density matrix takes the form 

p[$, $, T] = n A4(T) exp |-iA,(T) f^^{T) ■ f^_^,{T) - \aI{T) ~%{T) ■ ~ff^^,{T) 

k 

7r^(T) is the Fourier transform of n^(T, x). This form of the density matrix is dictated by 
the hermiticity condition 

p[$,l',T] =p*[<l,$,T] ; 

as a result of this, -Bfc(T) is real. The kernel -Bfc(T) determines the amount of 'mixing' in 
the density matrix since if = 0, the density matrix corresponds to a pure state because it 
factorizes into a wave functional depending only on $(■) times its complex conjugate taken 
at $(■). This is the case under consideration, since the initial conditions correspond to a 
pure state and under time evolution the density matrix remains that of a pure state |]13| . 



In conformal time quantization the evolution of the density matrix is via the conformal 
time Hamiltonian ( |6.5| ). The evolution equations for the covariances (given by equation 
(2.20) in see also equation (2.44) of [|^) are obtained from those given in ref. [|l3l by 
setting a(t) = 1 and using the frequencies uj^T) = k"^ + M.'^iT). In particular, by setting 



Kir) 

FkiX) ' 



10.2) 



More exphcitly 



MkiT) = Mk{%) exp 



r Aik{T') dT' 



A4(ro) 



AikiT) = -^\og\F,{r)\ 



d 



-a{t)-a{t) -\og\h{t)\ 



AukiT) 



1 



1 



|F,(T)P a{ty \Ut)\ 



2 ' 



Bk{r) = 



where Aji^ and Aj^ are respectively the real and imaginary parts of A^ and we have used 
the value of the Wronskian ( p.l2| ) in evaluating ( |10.3| ). 

The coefficients Ak{T) and Mk{T) in the gaussian density matrix (|10.1|) are completely 
determined by the conformal mode functions Fk{T) (or alternatively the comoving time 
mode functions fkif))- 

Let us study the time behavior of these coefficients. During inflation, a{t) ~ e^°^, and 
the mode functions factorize after horizon crossing, and superhorizon modes grow in cosmic 
time as in Eq. (p. 11): 



J2u~l)hot 



where the coefficient can be read from eq. (|9.1|) . 
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We emphasize that this is a result of the full evolution as displayed from the numerical so- 



lution in fig.|2l|. These mode functions encode all of the self-consistent and non-perturbative 
features of the dynamics. This should be contrasted with other studies in which typically 
free field modes in a background metric are used. 



Inserting this expression in eqs. (|10.3| ) yields 



AjkiT) -ho e'^«* ( - i ) + 0(e-^°*^ 



2 



Since z/ — | > 1, we see that the imaginary part of the covariance Ajk{T) grows very fast. 
Hence, the off-diagonal elements of $, T] oscillate wildly after a few e-folds of inflation. 
In particular their contribution to expectation values of operators will be washed out. That 
is, we quickly reach a classical regime where only the diagonal part of the density matrix is 
relevant: 

p[<l>, $, T] = n A4(T) exp {-AR,{r) v^,{T) r^_^,{r)} . (10.3) 

k 

The real part of the covariance Ajik{T) (as well as any non-zero mixing kernel Bk(T) 
||T3| ) decreases as Q-('^'^-'^)hot _ Therefore, characteristic fleld configurations r]^ are very large 
(of order e^'''~ ^^^°^) . Therefore configurations with field amplitudes up to 0{e^'^~^^'^'^^) will 
have a substantial probability of occurring and being represented in the density matrix. 

Notice that x ~ e^'^~2)'*o* corresponds to field configurations $ with amplitudes of order 



^{u ^)hot j-ggg eq. ( |6.2| )]. It is the fact that — | > which in this situation is responsible 



for the 'classicalization', which is seen to be a consequence of the spinodal growth of long- 
wavelength fiuctuations. 

The equal-time field correlator is given by 

(x(f,r) x(f',T)) = j |F,(r)|2 e^-^-'') , 

and is seen to be dominated by the superhorizon mode functions and to grow as ei'^'^-'^)^ot ^ 
whereas the field commutators remain fixed showing the emergence of a classical behavior. 
As a result we obtain 

(x(x,T) x(f' ,r)) oc a^{t) (peffit) 0e//(O G{\x-x'\) + small (10.5) 



where G{\x — x'\) falls off exponentially for distances larger than the horizon [T^ and 'small' 
refers to terms that are smaller in magnitude. This factorization of the correlation functions 
is another indication of classicality. 

Therefore, it is possible to describe the physics by using classical field theory. More 
precisely, one can use a classical statistical (or stochastic) field theory described by the 
functional probability distribution (|10.3|). 
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These results generalize the decoherence treatment given in ref. |^3| for a free massless 
field in pure quantum states to the case of interacting fields with broken symmetry. Note that 
the formal decoherence or classicalization in the density matrix appears after the modes with 
wave vector k become superhorizon sized i.e. when the factorization of the mode functions 
becomes effective. 



XI. CONCLUSIONS 

Since there are a number of articles in the literature treating the problem of preheating, 
it is useful to review the unique features of the present work. First, we have treated the 
problem dynamically, without using the effective potential (an equilibrium construct) to 
determine the evolution. Second, we have provided consistent non-perturbative calculations 
of the evolution to bring out some of the most relevant aspects of the late time behavior. In 
particular, we found that the quantum backreaction naturally inhibits catastrophic growth of 
fluctuations and provides a smooth transition to the late time regime in which the quantum 
fluctuations decay as the zero mode approaches its asymptotic state. Third, the dynamics 
studied obeys the constraint of covariant conservation of the energy momentum tensor. 

The Hartree approximation and the large limit are physically inequivalent. The 
Hartree approximation considers a single field with self-coupling g while in the large N 
limit one considers N ^ oo fields with an 0{N) invariant self-coupling g/N. Since grand 
unified theories contain a large number of scalar fields, the large N limit seems a more real- 
istic effective model than the Hartree approximation for the inflationary universe. Particle 
production during preheating turns to be much more effective in the large N limit than 



within the Hartree approximation (see sec. VII and [0]). In addition, the field evolution is 
definitely model dependent. For example, in the fine-tuned two-fields model (such that the 
Mathieu equation applies) studied in ref. parametric resonance in Hartree approxima- 



tion is significative only for a limited range of parameters. 

We have considered in these Lectures single-mass inflaton models. The Hartree approx- 
imation for a two-field model (with unequal masses) is presented in the second reference 
under [^. In such cases, one must carefully apply the Hartree approximation to all fields 
in the model self-consistently. Applying Hartree to one field and just doing perturbation 
theory in the other field leads to inconsistencies as noticed in ref. p[. 

It can be argued that the inflationary paradigm as currently understood is one of the 
greatest applications of quantum field theory. The imprint of quantum mechanics is every- 
where, from the dynamics of the inflaton, to the generation of metric perturbations, through 
to the reheating of the universe. It is clear then that we need to understand the quantum 
mechanics of inflation in as deep a manner as possible so as to be able to understand what 
we are actually testing via the CMBR temperature anisotropics, say. 

What we have found in our work is that the quantum mechanics of inflation is extremely 
subtle. We now understand that it involves both non-equilibrium as well as non-perturbative 
dynamics and that what you start from may not be what you wind up with at the end! 

In particular, we see now that the correct interpretation of the non-perturbative growth 
of quantum fluctuations via spinodal decomposition is that the background zero mode must 
be redefined through the process of zero mode reassembly that we have discovered. When 
this is done (and only when!) we can interpret inflation in terms of the usual slow- roll 
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approach with the now small quantum fluctuations around the redeflned zero mode driving 
the generation of metric perturbations. 

We have studied the non-equilibrium dynamics of a 'new inflation' scenario in a self- 
consistent, non-perturbative framework based on a large N expansion, including the dy- 
namics of the scale factor and backreaction of quantum fluctuations. Quantum fluctuations 
associated with superhorizon modes grow exponentially as a result of the spinodal insta- 
bilities and contribute to the energy momentum tensor in such a way as to end inflation 
consistently. 

Analytical and numerical estimates have been provided that establish the regime of 
validity of the classical approach. We flnd that these superhorizon modes re-assemble into 
an effective zero mode and unambiguously identify the composite fleld that can be used 
as an effective expectation value of the inflaton fleld whose classical dynamics drives the 
evolution of the scale factor. This identiflcation also provides the initial condition for this 
effective zero mode. 

A consistent criterion is provided to extract small fluctuations that will contribute to 
cosmological perturbations from large non-perturbative spinodal fluctuations. This is an 
important ingredient for a consistent calculation and interpretation of cosmological pertur- 
bations. This criterion requires that the model must provide many more than 60 e-folds 
to identify the 'small perturbations' that give rise to scalar metric (curvature) perturba- 
tions. We then use this criterion combined with the gauge invariant approach to obtain the 
dynamics of the Bardeen variable and the spectrum for scalar perturbations. 

We flnd that during the inflationary epoch, superhorizon modes of the Bardeen potential 
grow exponentially in time reflecting the spinodal instabilities. These long-wavelength insta- 
bilities are manifest in the spectrum of scalar density perturbations and result in an index 
that is less than one, i.e. a 'red' power spectrum, providing more power at long wavelength. 
We argue that this red spectrum is a robust feature of potentials that lead to spinodal in- 
stabilities in the region in fleld space associated with inflation and can be interpreted as an 
imprint of the phase transition on the cosmological background. Tensor perturbations on 
the other hand, are not modifled by these features, they have much smaller amplitude and 
a Harrison- Zeldovich spectrum. 

We made contact with the reconstruction program and validated the results for these 
type of models based on the slow-roll assumption, despite the fact that our study does not 
involve such an approximation and is non-perturbative. 

Finally we have studied the quantum to classical crossover and decoherence of quantum 
fluctuations by studying the full evolution of the density matrix, thus making contact with 
the concept of 'decoherence without decoherence' which is generalized to the interacting 
case. In the case under consideration decoherence and classicalization is a consequence of 
spinodal growth of superhorizon modes and the presence of a horizon. The phases of the 
mode functions freeze out and the amplitudes of the superhorizon modes grow exponen- 
tially during the inflationary stage, again as a result of long-wavelength instabilities. As 
a result fleld conflgurations with large amplitudes have non-vanishing probabilities to be 
represented in the dynamical density matrix. In the situation considered, the quantum to 
classical crossover of cosmological perturbations is directly related to the classicalization 
of superhorizon matter fleld modes that grow exponentially upon horizon crossing during 
inflation. The diagonal elements of the density matrix in the Schroedinger representation 
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can be interpreted as a classical distribution function, whereas the off-diagonal elements are 
strongly suppressed during inflation. 
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